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The subroutine naned VORTEX which is listed in the users manual 
has been modified to eliminate a defect which can cause errors when 
^win-engine aircraft configurations are modeled. 

Due to this modification* the following seven corrections must be 
made to the users manual. 

1. On page 21 in Table 2.1, the list of subroutine names in 
which Common Block /WlNGl/ appears must be changed. Remove the word 
VORTEX. Change that group of names to read 

MAIN PROGRAM, INPUT, COEFIC, 

COFSYM, WGFOM, WINGV, VELOCI, 
and VPROPS 

2. On page 25 under item d concerning Common Statement /WINGl/, 

delete the last line which reads: " -card VOR 2250 in subroutine VORTEX." 

3. On page 32 in line 17, change the number 3991 to 3975. 

4. On page 34 in line 8, change the number 363,782 to 363,766. 

5. On page 34 in line 26, change the number 25,266 to 25,250. 

6. In the Appendix A description of subroutine VORTEX found on 

pages 105 and 106, page 106 is correct. However, page 105 must be 
changed and should be replaced by the following attached page. 

7. Finally, the listing of subroutine VORTEX found on pages 225 
through 228 in Appendix B of the manual is no longer used. Instead, 
those four pages must be replaced with the listing of the newly modified 
subroutine VORTEX which is also enclosed on the last four pages of this 
errata. 

The sample case presented in Appendix C of the users manual is not 
affected by the subroutine modification and no corrections are needed 
there . 
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Savatnl caaaa of body angla of attack and aidcT allp nay ba aolvad 


.^ogathar during ona run of tha prograau Tha wing, if praaant, io 
"" f^praaantad by a alngla horsaahoa vortax bound at the ona qpartar 





ord lina. Tha trailing vorticaa ara aaparatad by a diatanea of 


ir/4 ciaaa tha apan. A vtog lift eoaffielant ia apaciflad by tha uaar 
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or aach angla of attack to tha wing vortax atrangth ia not an 


'^unknown in tha problan. Tha Influanca of tha wing on tha body ia 
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ponaidarad, but tha body la aaaunad not to influanca tha wit^. Tha 
body nay ba of any arbitrary ahapa. Tha body aurfaca it rapraaantad 
by a aaah of triangular and quadrilateral panala. Certain panala nay 
ba apaciflad at being ralaxad boundariaa through which tone apaciflad 
parcantaga of fraa atraam velocity ia apaciflad. Such panala nodal 
inlet and outlet ragiona. All other panala ara aolld boundariaa through 
fdiich no velocity nay paaa. 
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A cc'iatant diatrlbution of aourea atrangth acta on tha panel under 
conaidaration. At other panalt* tha aourea diatrlbution ia lunpad at 
a point aourea at tha panel control pointa. At tha panel under conald** 
aration* a aunnation ia nada of the noxnal velocity induced by the 
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1^ daalrad. foaat of eaxd output la ouch that tha carda aay ba uaad; 
diraetly aa part of tb* input data to c3m propallar parforaanca 
pradietion progran flvan in rafaraneaa 2 and 3. 
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This profraa ana writtan in aii^a praciaion FORXItAli for uaa on 
an IBM 370 conputar uaing an 08/370 oparatinf ayatan and uaing built- 
in library fynctiona* IWo auxiliary aaquai^ial aerateh filaa ara 
raquirad by tha progran. lhaaa filaa aara atorad on aa im 3330 Diafc 
Pack. Bataila on allocating apaea for thiwa filaa ia givan in faction 2 
of thia sanual. 
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A Hating of tha progran ia givan in i^^paadix B. Thia particular 
varaioa of tha progran ia dinanaionad to a lar^ alta eapabla of pro- 
caaaittg up to 2596 panala and 6 body oriantationt. Houavar, tha 
progran dinaaalona nay ba anlargad to anabla handling of vary long 
bodlaa, vary vlda bodiaa* or bodiaa raqulrlng nora panala or tha 
progran dlnanaiona can ba raducad to anabla tha progran to fit in a 
anallar conputar. Datailad procaduraa for altaring progran dintnaiona 
ara givan in Section 2 of thia nanual. Alao tha conqnitar atoraga 
raquirananta aa a function of progran dlnanaiona ia ai^lainad in Section 2 

Section 2 of thia nanual glvaa an in-dapth daacription of tha 
nachanicf of dila progran too detaHad to nantlon in thia brief 
daecrlption. 
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1«1 Pmi» 1 Dtstrlfcutleo a»i UtM Faeton 

Hm quality of and aueeasa obtaining raaulta from thla thraa* 
dlaanalonai potantlal flow prograa dapanda on tha nuaibar of body panala 
uaadt and tha dlatrlbutlon and alsa of panala. 

TWO factora ahould ba hapt In alnd «dwn paMllng tha configuration. 
The flrat factor, aa dlacuaaad In rafaianca 4, la concamad with tha 
dlatrlbutlon of panala. Tha panal aaah mat ba finer (nora and anallar 
panala) on raglona of tha body whara abrupt aurfaca curvature and rapid 
change In croaa aactlonal ahi^a occura. Exaaplaa of auch raglona are 
nacalla Inlet llpa, locatlona where falrlnga or other protruaiona are 
attached to tha body, and tha region where tha canopy or cabin inter- 
aacta tha fuaalaga. Kaarby gaomtry haa tha atrongaat Influence on tha 
propeller plana flow field. Since tha min purpoaa of thla prograa la 
to calculate tha flow field in tha propeller plana. It la daalrabla to 
noat accurately define tha raglona of tha body naaraat to tha propeller 
plana. Tharafora, apacial attao^ta ahould ba mda to accurately panal 
tha cowl or nacelle Inlet face and foraward body raglona uaing aa fine 
a paneling naah aa pcaalbla. Tha aacond factor to conaldar daala with 
panal aizaa. For beat raaulta, gradual chMgaa In panal alee nuat ba 
mda to bland panallng in ragima of danaa panal coneantratlona with 
panala in raglona of aparaa concentration. Rafaranca S atataa that tha 
charactarlatlc dlmnaiona of a panel ahould not differ by more than 
50 percent from tha dlmnaiona of adjacent panala. Although thla rule 
ia difficult to obey with every panal on all cmflgurationa. It glvaa 
a guideline for tha uaar to follow. 


Tha two factora of panal mah alze and dlatrlbutlon give tha uaar 
baaic guldaltoaa to follow in paneling tha body. Tha aourea atrangth 
aolutlona will generally converge to a aolutlon If tha two guidallnaa 
are followed. Tha accuracy of tha aolutlma obtained my depend aom- 
idtat upon how eloaaly tha guidallnaa arc followed. If tha body panallng 
atraya too far froa tha above guidalinea, tha mtrix ayataa of aquatlona 
my bacom ao poorly behaved (not diagonally dominant) that tha 
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lt«rativ« tolutloo proe*** will dlvtrgt giving no rooulcs. On eoopltx 
conf igurnciono , ospoeinlly mulciplt body typoo *• found in twin ongino 
Aircraft » tho paneling ntth nay load to dlvorgont aolutiona* In oono 
•ituationo* tho uoor miot uoo trial and error in aodifylng the paneling 
arrangenent to get convergent eolutione. Reference 5 contains a general 
diecueeion on the nuoMtical behavior of the eyaten of equatione charac- 
terietlc of thie paneling method. Other geometry paneling factors 
affecting solution convergence are explained in sections 1*2 and 1.3 
below. 

1.2 Sinaia Body Paneling Input Techniques 

1,2.1 Cross Section and Periphery Point Input Rules — General 

Exasiples of configurations considered as a single body are 
single engine aircraft cowl- fuselages without spinners, or the cowl- 
fuselage with the spinner attached to the cowling without a physical 
separation. Such a body is modeled by inputting the discrete point 
description of the body cross sections. See also section 3, card set 3 
in this manual. 

The method of paneling closely resembles the method used in refer- 
ence 4. The body is cut into many cross sections perpendicular or 
nearly perpendicular to the body reference x axis. Croes sections need 
not always be purely coplanar cuts through the body. Rather, the cross 
section cuts may \*reave or curve forward and aft in an Irregular fashion. 
See Figure 1.1. Such "quasi" cross sections may be needed to allow 
shifting of points forward and aft at awkwardly contoured regions on 
the body. 

Cross section descriptions mus^ be input In sequence starting at 
the front of the body end moving aft* See Figure 1.2. Note that a 
cross section may simply be a singl ,>oin'. repeated several times. 

This occurs in the example in Appendix C. An exception to the front 
to aft ordering occurs in paroling the inside of inlets (as at the front 
of a cowl). In this situatiot;, see Figure 1.3, the first cross section 
must be the one which defines the inlet contour and is farthest aft. 




Figure 1.2 Sample body cross section input showing a repeated cross section description 
and indicating the input sequence of periphery points. 
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The inlet wall cross sections are Input sequentially from aft to front 
to the Inlet lip. Then the outer contours are specified In the normal 
front to aft sequence. 

Each cross section Is defined by specifying the coordinates of the 
periphery points on the section in the sequential order Indicated In 
Figures 1.2 and 3.2. 

Between two cross sections, the program generates a ring of panels 
using the Input points for corners. See Figure 1.2. The ordering of 
the section points Is Important a? the first point on one section is 
paired with the first point in the next section, etc. By specifying 
a point twice or more, that point will become the comer of one or more 
triangular panels. Otherwise, quadrilateral panels will be generated. 

The input ordor of cross sections and periphery' points shown in 
Figures 1.2, 1.3, and 3.2 must be followed to ensure the panels will 
have outward dra\im normal vectors . 

Refer to Figure 1.2. Panel rings will be generated between cross 
sections I and I+l unless the flag signal, NFLAG, is given a value of 
1 on section I. NFLAG=1 for section I indicates that section I x^ill 
be repeated as section I+l using a different number or distribution of 
periphery points. Thus no panels will be generated between sections I 
and I+l in this situation. Use of NFLAG=1 is the method for increasing, 
decreasing, or altering the distribution of panels on the rings back 
along the body. It is essential that the number of points, NIP, on 
section I equal the same number of points on section I+l if panels are 
generated between these sections. 

Cross sections, between which panels will be placed, must never 
be allowed to intersect or coincide. That is, no point on the one 
section can be the same as a point on the other section. This would 
trigger an error in the program. However, cross sections may be con- 
centric. This might be necessary in paneling the flat front face of 
a cylinder, for example. 


It should be stated that although It is possible to panel concave 
body eotttouts, such as the Inside of an lnlet» It has been found that 
In auchcaaes, the equation solutions will usually diverge and no 
results will be obtained. It la not recoamended that concave ahapes 
such as the interior walls of inlets be paneled. Section 1.2.3 below 
presents an alternative technique for modeling Inlets. . 

1.2.2 Use of Symmetric and Non-Symmetrlc Input Options 

The Input of cross section periphery points completely 
around the section, as explained In section 1.2.1, is the general or 
non-synsetric Input option nethod (see Figure 3. 2.b). This input 
generates panels on both sides of the body and can always be used 
regardless of body symmetry. 

As is also outlined in section 3 of this manual, the symmetric 
body input option can be used If special conditions apply. These 
conditions are: 

1. The body cross sections must all have left-right symmetry 
about a comsion (x-z) plane of symmetry. This plane of 
symmetry must contain the y>0 coordinate. 

2. The root of the wing (if present) must be on this plane of 
symmetry. 

3. Inlets and outlets (panels) must be symmetrically arranged 
about the plane of symmetry. 

4. All body sideslip angles must be zero. The body may be 
oriented at angle of attack only. 

vnien all these conditions are met, then only points on the left 
(negative y) half of each cross section need be Input (see Figure 3.2.b). 
This Is the symmetric Input option. If this special Input option Is 
used, then the right side panels are mirror Images of those on the 
left and they are automat ally generated by the program. The 
loq>ortant benefit of the symmetric Input option Is that equations need 
only be created for the left side panels. The size of the system of 
equations Is half of that produced when non-symmetrlc Input Is employed. 
Thus computing time and costs are greatly reduced. It Is recommended 
that the symmetric Input option be used whenever possible. 
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Typically* a single engine aricraf t configuration can be paneled 
using the syinmetry option if the non-syametrical snail ducts and other 
small protrusions are Ignored* and if the spinner is attached to the 
fuselage symmetrically. If the spinner is yawed with respect to the fuselage 
(it typically is built this way)* then the symmetric input option can not 
be used with this configuration unless the spinner is not modeled. 

1.2.3 Relaxed (Inlet/Outlet) Panel Input Considerations 

. ■■ 

The body must always be paneled as a closed body. So if 
the front of the body is acttially an open inlet* the inside of the 
inlet must be paneled as a closed cup shape. However this paneling 
method is not recosmended as the solutions will fall to converge. 

Instead* the inlet face should be covered with "relaxed" panels which 
allow some flow to pass through them. To do this* panel the opening 
of the inlet as if It were closed off by a plate. Now in the inlet 
panel data cards (set 5 of input data— see section 3) the panel numbers 
of the inlet covering panels are specified. Assign a through flow 
velocity ratio to these panels. This velocity ratio will simtilate the 
flow into the mouth of the inlet. The configuration will appear as 
a convex closed body with inlet Inflow modeled. This should prevent 
solution divergence problems caused by concave paneled surfaces. 

If the user is uncertain of the inlet/outlet panel numbers to 
specify* he should first run the configuration using the geometry 
test run option* NCALC ■* 1. To do this* set NCALC * 1 on card set 4 
of input. Do not specify any inlet/outlet panels in card set 5 of 
input. Submit the data. The program will generate the body paneling 
but will not solve the flow equations. Take the output panel geometry 
from the run and find the numbers of the panels which should be inlets 
or outlets. Go back to the input data and specify these inlet/outlet 
panels in card set 5. Change the NCALC > 0 on card set 4. Then the 
data is resubmitted and a normal run will be made. All flow calcula- 
tions will be made. Refer also to section 3 of this manual for use 
of the test run* NCALC > 1* option. 
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eoafiiuration nay consist of several discrete 
^osed iiodles. An exao^le of this la a single engine aircraft fuselage- 
spinner edgblnatlon In which the spinner Is paneled as one Isolated 
closed body and the fuselage Is a second Isolated closed body. A second 
esanple treated as a multiple body Is e twin engine aircraft configura- 
tion with wing mounted nacelles. Each nacelle Is a separate closed 
body. The fuselage Itself Is a third closed body and Is not physically 
connected to the nacelles. 

Multiple body configurations are paneled by considering one body 
at a time using the techniques given In section 1.2.1. All bodies 
are referenced to one common coordinate system. It does not matter 
which discrete body Is defined first. For exaople» a nacelle can be 
paneled first, then the fuselage, and lastly another nacelle. Or the 
fuselage could be paneled first. If, on a twin nacelle aircraft, the 
left propeller plane Is considered for the predictions. It may be 
convenient to input the left nacelle body first. 

Input the first discrete body cross sections according to 
section 1.2.1 rules. Use the NFLA6 * 1 parameter If a cross section 
repeat is desired. To distinguish the end of the first discrete body 
from the start of the second discrete body, It is necessary that 
MFLAG ■ 1 be specified on the last Input section of the first body. 

In this use of NFLAG ■ 1, the last cross section of the first body 
will nor be repeated. Instead, the flag signifies that no panels will 



b« placed between thla leet eroee eectlon of the firat bodf end the 
flret croee eectlon of the eecond body. Now Input the eecond discrete 
body using the rules of section 1.2.1. \s before. If e third body 
follows the second, use NFLA6 • 1 on the lest cross section of the 
second body. Panel the third discrete h^dy etc. Panel the lest 
discrete body using the rules of section 1.2.1. However, since the 
last discrete body will have none following It, do not specify NFLAG*1 
on the last section of the last body. The result Is several closed 
bodies, each paneled according to single body rules with none of the 
bodies physically connected. 

Through experience In paneling nacelle>fuselage-nacelle configura- 
tions, It has been found that this panel geooetry often produces a 
system of equations which will diverge when solved Iteratively. To 
obtain results. It may be necessary to Ignore some regions of the con- 
figuration. For example. If only the left nacelle and foreward portion 
of the fuselage are paneled, the equations may be well behaved and will 
converge. Then if flow calculations are made on the left propeller 
plane, the results are still valid since the left nacelle and front 
fuselage have a dominant influence on this propeller plane. The right 
nacelle and rear fuselage are remote from the left propeller plane and 
it is safe to Ignore these components in paneling the configuration. 

1-3.2 Use of Symmetric and Non-Symmetric Input Options 

Each discrete body of a multi-body configuration may always 
be input using the non-symmetric input option. To use this option, 
follow the procedures in sections 1.3.1 and 1.2.1 (see Figures 1.2 and 
3.2. a). 

The same conditions as given in section 1.2.2 must hold before the 
symmetric body input option may be used on multi-body configurations. 
These conditions are repeated here and have new meaning when applied 
to multiple body configurations: 

1. All cross sections of all discrete bodies must be symmetric 
(left and right) about one common plane of symmetry. This 
x-z plane of symmetry must contain the y ■ 0 coordinate. 


2. Th* root ceordlnat* of tho wing (If prooont) gw«t llo on this 
^ pisns/of 'synastry. ' 

Inlsts and outlsu (pansls) oust bo syiaBstrlenlly srrangsd ! 
id^ttt, ths ; plans ;(^,,j#yai»try< 

4. All body aldssllp angles must b« aero. The configuration 
tti^ be orientad at angle of attack only. . ^3 ; 

All of these four conditions iDust exist before the synastric 
input option can be used. In using this syonetric option, only the 
points on the left (negetive y) side of each body cross section ere 
input (see Figure 3.2.b). The right side nirror inage panels will be 
autooetlcally created by the prograa. However, equations need only 
be written for the left side panels. Thus, the system of equations is 
only half as large as would occur if the non-synmetric option were used. 
Savings in cosputer time and cost is the result. 

Note that rule 1 above requires that all discrete bodies of a 
multiple body configuration oust lie in tandem on the common x axis. 

This is the only manner by which all cross sections of all bodies could 
be symmetric about one plane of sysMtry. A multiple body configuration 
which could use the symmetric input option is a single engine fuselage^ 
spinner coid>ination, provided the spinner is not yawed with respect to 
the fuselage. 

Although typical twin engine nacelle-fuselage-nacelle configura- 
tions would seem to be symmetric, they can never be paneled tislng the 
symmetry option as defined for this program, because each nacelle lies 
completely to one side of the plane of symmetry and does not obey rule 1 

Other examples of configurations which can not use the symmetry 
option are: 

- nacelle-nacelle bodies side by side (rule 1 does not hold) 

- fuselage-yawed spinner (rule 1 does not hold) 

- fuselage-nacelle (rule 1 does not hold) 

- wing mounted nacelle and wing. This is a single body which 
cannot use the symmetry option because the wing root does not 
lie on the nacelle plane of symmetry. 




Uoforeunatcly, at aaan abova, eha tynBMtrie input option haa 
llnittd application whan panallng aulclpla body configurationa. 

1.3.3 ftalaxad (Inlat /Outlet) Patial Input Cenaidarationr, 

Each diacrata body of tha suit ipla body configuration auat 
ba paoalad aa a cloaad body. If tha bodiaa contain inlat or outlat 
opaninga, tha uta of ralaxad inlat/outlat panala can ba sada. Thaaa 
ralaxad panala can ba uaed to sodal inlat opaninga without panallng 
tha actual concave innar wall ahapa. Tha dlacuaaion ox thaaa ralaxad 
panela given in aactlon 1.2.3 appllaa equally to aultlple body 
configurationa. 

1.4 Sunmary 

The forsulaa uaad for calctilatlng panel charactarlatica fros 
diacrata croaa aactlon periphery point input can be found in referanco 1. 
This aection of tha manual ahould be used in conjunction with aactlon 3 
when Baking the input data deck. 

Although any arbitrary body geometry can ba paneled cucceaafully* 
the syataa of equationa obtained from the paneling network may not be 
well behaved numerically. As a result, possible solution divergence 
nay occur during the iterative solution process. Thus, it may be 
necessary to use trial and error in getting a paneling network which 
accurately represents the body and also converges to a solution. 
Basically, the quality of the paneling depends upon the Judgement and 
experience of the user. 



SECTION 2 OVERALL PROGRAM DETAILS 

;w 

2*1 Proaran Structurt 

The 3-D potential flow cooputer program conaiata of a main program 
uaed in aettlng up array dlmenalona, eatabliahing program alae, and 
calling the aubroutines. Of the 14 aubroutlnea uaed In the program, 

10 are called directly by the t?aln program. The overall atructure of 
the program la illuatrated by the flow charg on Figure 2.1. Alao a 
dcacrlption of the main program and each aubroutlne la given In 
Appendix A. 

A copy of the program la Hated In Appendix B. Thla particular 
veraion haa been dimensioned to a fairly large size capable of handling 
up to 2596 body panels and up to six body orientations. The size of 
the program In Appendix B la one of the largest which can be compiled 
and executed on the IBM 370 computer facility at the Pennsylvania State 
University. 

2.2 Program Dimension Size — Introduction 

As the program dimension size determines Che size of the configura- 
tion which may be handled, core storage and auxiliary file storage 
requirements. It may be necessary to change the size of the program 
given In Appendix B. Some or all of the dimensions of Che program may 
be decreased to enable use of the program in a smaller cosq)uter Installs 
tlon, or to enable the program to be run In a smaller storage-higher 
priority run category at the Installation. Conversely, some or all 
dimensions of the program may be Increased to take advantage of a larger 
computer facility. This would enable configurations with more panels to 
be handled. 

In Che remaining portions of this section is a description of Che 
variables which must be changed by the user to alter the program 
dimensions. Also described are the C0M^K}N statements and DIMENSION 
statements needed, auxiliary files, and step-by-step procedures to be 
followed in altering the dimension size of the program. Then instruc- 
tions for allocation of auxiliary file storage are given along with 
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ectlmations of output record rcquireaents end core storege requlronence 
for the given progran difflenelone. Execution time le dlecuaeed, and 
Inetructlona are given for creating an object program deck to enable 
the program to execute faster. Lastljr* a section gives Instructions 
for operating the program, Includii^g all required j b control cards 
needed to run the program on the Penn State compute facility. 


2.3 Program Dimension Size Controlling Variable Descriptions 


The user has control of the progran dimensions by specifying 
values of the following four variables on four cards In the main 
program. The main purpose of these is to permit execution time 
dimensioning of certain arrays In the subroutines: 


1. NSECMA Upper i’imlt on the number cf body cross section 

descriptions which may be input to define the body. 
Rule : NSECMA must be integer value of 3 or more. 

2. NIPMAX Upper limit on Che number of periphery points allowed 

around a con^lete cross section. (First and last 
point are the same but count as two points.) 

Rule ; NIPMAX must be an odd Integer with a value of 
S or more. 


3. MAXALP Maximum number of body orientations (pairs of body 

angle of attack and sideslip) which may be Input in 
one run of the program. 

Rule ; MAXALP must be an integer value of 1 or more. 

4. MAXINF Maximum number of body panels over tii? entire configura- 

tion which may be specified as being inlet or outlet 
panels. 

Rule ; MAXINF must be an even integer value of 2 or 
more. 


lx the user breaks any of the rules above, the program will detect 
them and stop. Refer to section 2.3 for instructions on changing these 
values, and refer to the listing in Appendix B as an example. 


18 


« 

Ba««d on the u««r apaelflad «lza control varlablaa akovt, chraa 
other variables are eooputed by the prograa and are also used for 
execution tiae diaenslonlng of certain arrays in subroutines. However, 
the user must also manually calculate the values as he must supply the 
nu^rs in the oain program DIMENSION statesMnt (saction 2.5) and in 
the COMMON statements (section 2.6). The three additional variables 
are: 

5. LLL LLL • (NSECMA - 1)*(NZPMAX - 1). LLL is the maximum 

possible nuo^er of body panels generated. This amount 
can occur only if NSECMA cross sections have been input 
with NIPMAX points on all sections. Also the paneling 
input would have to be completely efficient, that is, 
no repeated deserlptlotti of a cross sections and no 
triangular panel input. Cenerally, with complex 
bodies, cross sections are often repeated with two 
descriptions and section points are repeated several 
times to create triangular panels. So LLL panels will 
uot usually be generated, but could be. 

Note ; LLL will always be an even Integer of size 8 
or mere. 

6. MD MD ■ (LLL MAXALP) . MD is the maximum number of 

terms (right and left side coefficients) in each equa* 
tlon for panel source strengths. MD also is Che length 
of logical records which must be specified in auxiliary 
file 9 (see sections 2.7 and 2.9). 

7. LLLKAF LLLHAF ■ (LLL/2). It is the maximum possible number 

of panels on the left half of a symmetrically input 
body. 

See section 2.8 /.or the procedures involving variables LLL, MD, 
and LLLUAP . 
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2.4 DPSlp^^^ anU NRADj^j^^—Ptflnltiong 

iVo input dact valuta control the nuabtr and diatributlon of flow 
aurvty polnta in the proptlltr plant. Tht aininun value of one and 
toaxinun value of the other control the aize of auxiliary file 50 and 
the aaount of output recorda: 

1. DPSl AzitButhal angular increment between radial apokea of 

aurvty polnta on the propeller plane. DPSI b^^ ia t’^e 
saalleat increment allowed and nuat be more than 0.0 
but not more than 360. A value of 1.0 ia a 

good choice and ia u»ed in the program veraion in 
Appendix B. 

2. NRAO Nuii^er of aurvey polnta poaitioned radially from hub to 

tip at each azimuth increment on the propeller plane. 

It muat be a poaitlve integer but not more than NRABqhx* 

An NRAD^^ value of il la a good choice and haa been 

uaed in the progrra of Appendix h. 

>> 

The choicea of DPSIoin determine the maximum number 

of aurvey polnta generated. Thua they affect the number of output 
recorda produced. More importantly, they determine the I'.uii^er of 
records which muat be provided in auxiliary file 50 (aee aectiona 2.7 
and 2.9). 

The user must know what theae valuea are in the program version 
being used sc he cen properly allocate file 50 storage. Instructions 
for changing DPSIair. '-^^max values in the program arc found in 
section 2.3. 

2.5 Description of Array DIMENSION Statement In Main Program 

The following array names are found in a DIMENSION statement at 
the beginning of the main program. Refer to Appendix B. The dimen- 
sions of these arrays must be changed by the user in this statement 
when the program dimension size is altered (see section 2.8, step 3). 
These arrays are the only ones which are automatically dimensioned 


within the subroutines during execution. The DIMENSION statement end 
arrays are given below (without subscripts): 

DIMENSION A, SIGMA, NCOUNT, NFLAG, PX, PY, P2, STOXl, STOYl, 

STOZl, STOX2, STOY2, ST0Z2, ANVX, ANVY, ANV2, S, XC, YC, ZC, 

BIG, VARIAB,- SUM, SIGSAV 

2.6 COMMON Block Descriptions 

A total of eight different COMMON blocks are used throughout the 
program. Some subroutines use most of the COMMON statements and others 
use but a few. In most cases these COMMON statements contain both 
subscripted (array) and unsubscrlpted variables. The array variable 
dimensions In the COMMON statements must be changed by the user when 
the program dimension sizes have been changed. Dimension changes must 
be made to every COMMON statement card found in the program. The 
procedure for making the changes is found in section 2.8, step 4. 

Table 2.1 describes each COMMON block name and indicates the type 
of data stored and the locations at which the COMMON block cards must 
be present. 

2.7 Auxiliary File Descriptions 

File 3 is a formatted card data input file. File 6 is a formatted 
output file written on a standard 132 character-per-line printer. 

File 7 refers to a formatted data output file punched on standard 80 
column computer cards. These three file numbers are the default values 
for the computer facility at the Pennsylvania State University. 

Additionally, tx^o auxiliary scratch files are needed by the 
program. They are files 9 and 50 and nay be stored on magnetic tape 
or disk. These files are described below: 

File 9: This is a sequential scratch file used only to store the left 

and right side coefficients of the system of linear equations 
for the unknown source strengths. Each logical record of 
this file contains the coefficients of one equation. The 
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records save the equations in sequence starting with the 
first equation or row of the system. 

This file is written during execution of subroutine 
COEFIC if the symmetric body input option is not used. Then 
one record (equation) is written for each body panel. If the 
symmetric body input option (see section 3) is used, the 
file is written during execution of subroutine COFSYU. In 
this case, a record (equation) is written only for the panels 
on the left side of the symmetric body. 

The file is read during the Iterative solution for panel 
source strengths during execution of either subroutine SOLVE 
or SOLSYM. 

On the Penn State computer facility this file is stored 
on an IBM 3330 Disk Pack. 

File 50: This is also a sequential scratch file which temporarily 

saves propeller plane flow field prediction output. This 
allows the output to be sequentially printed on propeller 
plane output t.able part 2 (see section 4). The file is 
written and read only during execution of subroutine VPROPS. 

Each logical record of file 50 contains flow output 
information for one survey point in the propeller plane. 

Each record of the file corresponds to one line of the printed 
output table. Data written on each record consists of nine 
values including radial and azimuthal location of the point, 
three velocity components and four flow angles. 

On the Penn State computer facility, file 50 is scored 
on an IBM 3330 Disk Pack. 

The details of allocating file 9 and file 50 storage space is 
given in section 2.9. Also, that section ejqjlains how to write the 
job control cards for these files for use on the Penn State facility. 
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I 2.8 User Procedure for Altering Program Dimension Size 

If the user wishes Co alter the dimensions of the program from 
chose of Che version listed in Appendix the following steps must be 
performed. In what follows refer to sections 2.3 to 2.6 for the defini'- 
tions of dimensioning variables » dimensioned arrays and C0Mifi3N blocks. 
Refer to the program listing In Appendix B when a specific card number 
is mentioned in Che steps. Here are the nine steps involved in 
modifying Che program dimensions: 

Step 1. Choose the desired values of dimensioning variables 

NSECMA, NIPMAX, MAXALP. and MAXINF. Ensure they conform 
to the rules given for them in section 2.3. 

Step 2. Before changing the program cards, check the new core 

storage requirements of the program corresponding to the 
values chosen in Step 1. Use the procedure outlined in 
section 2.11 for this. If it is found Che program will 
require more storage Chan is available, it will be 
necessary to repeat Step 1 with new values until Che 
program core storage required is acceptable. Then perform 
the modifications given in the remaining steps 3 to 9. 

Step 3. Physically replace cards MAN 5850, MAN 5900, MAN 5950, 
and MAN 6000 of the main program with cards having the 
new values of NSECMA, NIPMAX, MAXALP, and MAXINF. 

Step 4. Calculate the values of the other three sizing variables 
LLL, MD, and LLLHAF using the values of Step 1 in the 
following formulas: 

a. LLL * (NSECMA - 1)*(NIPMAX - 1) 

b. >© - LLL + MAXALP 

c. LLLHAF - (LLL/2) 

These values are automatically generated by the program, 
but the user must know the values so the DIMENSION state- 
ment in the main program can be modified in the next step. 
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Ifodlfy the DIMENSION stateiDent in the main program by 
physically changing cards MAN 5200 to MAN 5350. Qiange 
the array dimensions to values compatible with the chosen 
dimension sizes of Steps 1 and 4. To do thls» use the 
following pattern DIMENSION statement which shows the 
array dimensions arbitrarily in terms of the dimensioning 
variables. Repunch cards MAN 5200 to MAM 5350 using the 
pattern below. (Insert the actual number values of the 
array subscripts, however. Punch the characters in the 
columns indicated. Also punch the card identifier numbers 
In columns 73 to 80.) (Note: only four cards should be 

needed for this statement.) 

Card Column: 

789 72 

DIMENSION A(MD) ,SIGMA(LLL,MAXALP) ,NCOUNT(N3ECMA) , 
NFLAG(NSECMA) ,PX(NSECMA,NIPMAX) ,PY(NSECMA,NIPMAX) , 
PZ(NSECMA.NIPMAX) ,ST0X1(NIPMAX) ,ST0Y1(NIPMAX) , STOZl (NIPMAX) , 
STOX2 (NIPMAX) , ST0Y2 (NIPMAX) , ST0Z2 (NIPMAX) , ANVX(LLL) , 
ANVY(LLL) ,ANVZ(LLL) ,S(LLL) ,XC(LLL) ,YC(LLL) , ZC(LLL) , 
BIG(MAXALP) ,VARIAB(MAXALP) .SUM(MAXALP) ,SIGSAV(ELLHAF) 

Five of the COMMON blocks contain subscripted array 
variables which may be changed. Modify the dimensions 
of these COMMON statement arrays affected by the changed 
dimension size variables of Steps 1 and 4. Every affected 
COMMON card in the main program and subroutines must be 
changed by using the new numerical values of the dimension 
size variables. The COMMON blocks which must be changed 
are given below and the cards to be changed are also 
given. (Note: where the dimension size variable names 

appear in parentheses, they are to be replaced with the 
actual numerical value). Also the COMMON statements 
must be pu.iched with the variables exactly in the sequence 
shown: 
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a, COMMON /INLET/ FLRATO(MAXINF) ,INDEX(MAXINF) .NINFLO 


-card INP25450 
-card COF 1300 
-card CSY 900 
-card VEL 1150 
-card VFR 2450 


in subrojtlne INPUT 
In subrov'Clne COEFIC 
In subroutine COFSYM 
in subroutine VELOCI 
in subroutine VPROPS 


b. 


COMMON /INPUTS/ ALPHA(ilAXALP) .BETA(MAXALP) , 
VX(MAXALP) ,VY(MAXALP) .VZ(MAXALP) , V.NALPHA 


-card MAN 5400 in 
-card INP25200 in 
-card VCM 1100 in 
-card COF 1250 in 
-card SOL 1200 in 
-card CSY 850 in 
-card SOS 1350 in 
-card VEL 1100 in 
-card VPR 2200 in 
-card WGE 2400 in 
-card VOR 2200 in 


main program 
subroutine INPUT 
subroutine VCOMP 
subroutine COEFIC 
subroutine SOLVE 
subroutine COFSYM 
subroutine SOLSYM 
subroutine VELOCI 
subroutine VPROPS 
subroutine WGEOM 
subroutine VORTEX 


c. COMMON /SYMTRE/ NPNSYM(LLL) ,INSOLV(LLLMF) ,NSYMET 

-card MAN 5450 in main program 
-card INP25400 in subroutine INPUT 
-card PAN 1300 in subroutine PANEL 
-card CSY 950 in subroutine COFSYM 
-card SOS 1400 in subroutine SCLSYM 


d. COMMON /WINGl/ CL (MAXALP) , CHORD, DIKED, SPAN, SWEEP, 
XQR,YQR,2QR,NWIMG 


-card MAN 5500 in 
-card INP25500 in 
-card COF 1350 in 
-card CSY 1000 in 
-card WGE 2450 in 
-card VEL 1200 in 
-card VPR 2400 in 
-card WGV 2800 in 
-card VOR 2250 in 


main program 
subroutine INPUT 
subroutine COEFIC 
subroutine COFSYM 
subroutine WGEOM 
subroutine VELOCI 
subroutine VPROPS 
subroutine WINGV 
subroutine VORTEX 


e. COMMON /WING2/ GAMMA(MAXALP) ,XTRALL(MAXALP) ,YTRALL 

(MAXALP) ,2TRALL (MAXALP) ,XTRALR(MA:iALP) ,YTRALR(MAXALP) , 
ZTRALR(MAXALP),XBTIPL,YBTIPL,ZBTIPL,XBTIPR,YBTIPR,ZBTIPR 

-card WGE 2500 and WGE 2550 in subroutine WGEOM 
-card WGV 2850 and WGV 2900 in subroutine WINGV 


! 


Step 7. Modify the following cards in the main program and sub- 
routine INPUT as explained below: 


a. 


b. 

c. 

d. 


Main Program 

card MAN 4400: insert new values of NSECMAyNIPMAX, 

MAmP, and MAXINF 


card >1AN 4550: Insert new value of LLL after the 

words "UP TO LLL 

card MAN 4700: insert value of MAXALP after the 

words "UP TO" 


card MAN 4800: insert value of MAXINF after the 

words "UP TO" 


e. 

f. 

h. 

i. 


j. 

k. 

l . 


m. 

n. 


Subroutine INPUT 


card INF 1850: 
card INP 2300: 
card INP 4050: 

card INP 5650: 
card INP 6650; 

card INP 6900: 
card INP11200: 

card INP11900: 

card INP12000: 
card INP14050: 


"C . . . CONTAINS 3 TO N CARDS ..." 
insert value of N where N ■ (MAXALP+2) 

"C . . . (VALUE FROM 1 to N . . 

Insert value of N, where N » MAXALP 

"C . . . CONTAINS 11 TO N CARDS ..." 
insert value of N, where 
N » (NIPMAX * NSECMA) + NSECMA + 2 

"C . . . MAX OF N) . . ." 
insert value of N, where N NSECMA 

"C . . . N). IF . . ." 

insert value of N, where 
N = (NIPMAX + l)/2 

"C . . . MAXWtUM OF N . . ." 
insert value of N, where N => NIPMAX 

"C . . . CONTAINS 1 TO N CARDS) ..." 
insert value of N, where 
N » (MAXINF + 1) 

"C . . . VALUE N . . 

Insert value of N, where 
N - (MAXINF/ 2) 

"C . . . VALUE OF N IF . . ." 
insert value of N, where N » MAXINF 

"C . . . N CARDS) ..." 

insert value of N, where N » MAXALP + 2 
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Step 8. If desired, change either or both the values of DPSl^ln 
and NRADg)ji 2 (. Ensure they obey the rules for these values 
given In section 2.4. The reconmended values of OPSIoin <■ 
1.0 degree and NRAD ^av ** 51 have been used In the program 
listed In Appendix B. If these values are changed, perform 
the following: 

a. If DPSImin Is changed, modify the following card 
found In subroutine INPUT 

-card INP20300: "C . . . MUST BE AT LEAST N AND . . .” 

Insert value of N, where 
N » DPSIjjin value 

b. If NRADjjax Is changed, modify the two following cards 
found In subroutine INPUT: 

-card INP21100: "C . . . FROM 1 TO N . . ." 

Insert value of N, where 
N » NRAD „,ay value 

-card INP36450: ". . .20 IF (NR.iJ).GT.N) GO TO 76" 

Insert value of N, where 
N - NRAD^ax value 

Step 9. The last step In modifying program dimensions is this. 

Calculate the amount of auxiliary storage needed for 
files 9 and 50. Write the job control language cards to 
give the new storage allocation. The rules for this are 
given in section 2.9. 

The version of the program given in Appendix B has been dimensioned 
using the following dimension size variable values: 

NSECMA - 60 

NIPMAX • 45 

MAXALP » 6 

MAXINF » 500 

LLL » 2596 
I-fl) » 2602 

LLLHAF « 1298 

and 

NRADtnax * 51 

DPSIroin = 1.0 
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2 . 9 Determining Auxiliary File Storage. Allocation Based on Program 

Dimenelon Size 

Storage for auxiliary file 9 and file 50 must be allocated on 
either magnetic tape or disk, and job control cards must be made to 
perform the allocation. The maximum dimension of the program version 
being used deterTnines how many records must be available on the two 
files. This was explained in section 2.7. 

In the following parts 1 and 2 below is given the formulas for 
calculating record and disk pack requirements for files 9 and 50, 
respectively. Also given is the job control language needed for 
establishing these files on the Penn State computer facility. 

2.9.1 File 9 Storage Allocation Formulas and JCL 

File 9 is a sequential fils requiring variably spanned 
physical records^. Enough cape or disk space must be provided to 
store up to LLL (section 2.3) logical records, each containing 
MD (section 2.3) single precision numbers of 4 bytes each. 

On the Penn State University IBM 370 computer system, file 9 is 
stored on an IBM 3330 Disk Pack^. Since this file is relatively 
large, Che disk storage is allotted by requesting numbers of cylinders 
of space, CYL. 

For this disk pack, the user should specify the following param- 
eters in the job control cards when allocating space for file 9; 

RECFM « VS (i.e., variable spanned records) 

BLKSIZE » 3120 (bytes) 

LRECL - (BLKSIZE - 4) - 3116 (byt?s) 

CYL ■ number of disk cylinders requested (calculated below) 

Number of cylinders, CYL, required is a function of the size 
variables, LLL and MD. Use the formula 1 or 2 below (whichever 
applies) : 
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1. If [(4 * HD) > (BLKSrZE - 8)] use; 

MD * 4 


CYL 


LLL * 


19 * 




BLKSIZE •> 8^ rounded to 
higher Integer 


1 

J 


13165 

BLKSIZE * 135 


rounded to 
lower integer 


^rounded to 
higher Integer 


-or- 

2. If ((4 * >D) :• (BLKSIZE - 8)) use: 


CYL 


19 * (: 


LLL 

13165 


t). 


BLKSIZE + laS-* rounded to 

lower Integer 


j rounded to 
higher Integer 

Note that CYL Is the niaximum number of cylinders which would be 
needed If LLL body panels were generated. So CYL must be available, 
but usually, not all of the space Is used, since most configurations 
seldom generate the full amount of LLL panels. 


The following example specifies file 9 storage and job control 
cards for the program version listed in Appendix B. This example 
applied to the Penn State University computer facility. This program 
has: 

LLL - 2596 and - 2602 
The job control parameters are; 


BLKSIZE - 3120 

LRECL - 3116 

CYL » 137 calculated using formula 1 above. 

Two job control cards punched exactly as shown below, both starting 

in column 1, will establish 137 cylinders of space for file 9; 

//DATA.FT09F001 DD UNIT-SYSDA,SPACE-(CYL. (123,1) .RLSE) . 

/ / DCB» ( RECFM-VS , LRECL- 3116 , BLKS IZE- 3120 ) 
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Here Che cylinders have been requested using a prlaary allocation 
of 123 cylinders with a secondary allocation of one cylinder (repeatable 
up to 14 times) which Is made when the space Is needed. This gives 137 
cylinders available. Reference 6 contains details on Che disk pack and 
Job control language for the Penn State computer system. 

2.9.2 File 50 Storage Allocation Formulas and JCL 

File 50 Is a sequential file requiring variably spanned 
physical records^. Enough cape or disk space must be provided to score 
up to R logical records, where R Is given by: 

» ' ‘ (Bii^)round.d to 

lower integer 

NRAD and DPSI , are defined in section 2.4. 
max min 

Each record contains nine single precision numbers of 4 bytes each. 

On the Penn State University IBM 370 .computer system, file 50 is 
stored on an IBM 3330 Disk Pack^ . The disk storage for this file 
should be allotted by requesting numbers of cylinders of space, CYL. 

For this disk pack, the user should specify the following 
parameters In the Job control cards when allocating space for file 50: 

RECFM ■ VS (i.e., variable spanned records) 

BLKSI2E - 44 (bytes) 

LRECL ■ (BLKSI2E - 4) • 40 (bytes) 

CYL ■ number of disk cylinders requested (calculated below. 

Mumber of cylinders, CYL, required is a function of tv.'^ propeller 
plane input point limits, NRi\D „^a y and OPSI^jin# of the program version 
being used. Usi the following formula to calculate CYL for file 30: 
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CYL 


NRAJWx (■ 


360 




DPSlg^n^ roundtd to 

lower Inf gT 


19 * (l 


13165 


r). 


BLKSIZE + 135-' rounded to 

lower integer 


rounded to 
higher Integer 


The following exaoq)le specifies file 50 storage allocation .ind Job 
control cards for the program version listed In Appendix B. This 
example appllea to the Penn State University computer facility. This 
program has the following propeller plane survey point limits: 

DPSI^n ■ 1*0 NRAD .mv ■ 51 


The job control pa;;ametera are: 

BLKSIZE > 4A 

LRECL - 40 

CYL ■ 14 calculated using the above formula. 

TWO job control cards punched exactly as shown below, both starting 

in column 1, will establish 14 cylinders of space for file 50: 

//DATA,FT50F001 DD UNIT-SYSDA, SPACE- (CYL, (1,1),RLSE), 

// DCB-(RECF:>VS,LRECL-40,BLKSIZE-44) 

Here the 14 cylinders have been requested using a primary alloca- 
tion of one cylinder with a secondary allocation of one cylinder 
(repeatable up to 14 times) which Is made when the space Is needed. 

This gives 15 cylinders available. Reference 6 contains details on the 
disk pack and job control language for the Penn State computer system. 


2.10 Output Record Considerations baaed on Program Dimension Size 

For some computer installations such as the one at Penn State, 
there Is a limit on the number of output records which may be produced 
during a single run. At Penn State, this output Limit Is 30,000 records 
If more records will be produced, special handling must be made to have 
the records written to a tape so It can be printed later In smaller 
portions . 
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Th« ntxlBun nunbtr of output rocords froa thij coisputor prograa io 
a ftnction of th« dlaaniion aiaa varlablaa (aaction 2.3) and ptopallar 
plana Input Units (section 2.4). The four primary quantities 
affecting length of output are LLL, MAJUl.?, 

An approxinate formula for estimating the aaxiiMin nusiber of output 
records produced is given below. The formula is a function of the 
dimension size variables defined earlier in sections 2.2 and 2.4. Also 
used are certain input variables defined in section 3. Note the formula 
assumes printed output is written on a standard 132 character-per-line 
printer* and the punched output is on SO column cards. Refer also to 
section 4: 

i- R^ + Rj + Rj Rj + MAXALP * (R^ + R^ + Rg + R^) 


where R^^^ is the total maximxin possible number of output records 
produced by a given size program version. 


The R^ terms in the formula are each given by formulas below: 


Ri 


0 * if. program listing is not printed 

3991, if program listing is printed * (source ..istlng) 


MAXALP + 1, if NWING « 1 
0, if SWING * 0 



if NPOINT - 1“| 
if NPOINT • 0 


■►MAXALP 


■► MAXINF NSECMA * (NIPMAX -H) ■»• 27 : 


(This is list of input.) 


R 


3 


0, if SWING « 0 

(MA.XALP •►26), if SWING • 1 


(This is wing geometry cable.) 


^ , T., if SLIST • I and NCALC • 0 
^ (5 * LLL) •►16, if NLIST - 0 or NCILC • 1 


(This Is body 
geometry table.) 


Re 


0, if yCALC - 1 

(3 * UMAX) •► 4, if NCALC - 0 


(This is solution iteration 
table. ) 


^ ^ jO, if sr.\LC • 1 

^ I (LLL + 14), if NC\LC - 0 
L 


: (Surface velocity table.) 
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■ 4 - 


•« 


0, It HCALC • 1 

mms, it N?0INT - 0 Md ilCAlC • 0 

«360/0PSI,la) . . * HllAB^l ♦ M, if NPOINT . 1 

lnM,« «d UCAIX - 0 

(ihls !• propeller plane oucpuc table part 1.) 


fo. If 

UiW.t, 


MCALC - 1 


R. - HmiP, if NPOINT - 0 and NCALC - 0 


|U360/DPSIain)rounded to * NRAD^xl + 3^. if J«WINT - 1 
" lower integer ‘ ° 

(Thie ie propeller plane output table part 2.) 

[b, if NPUNCH ■ 0 

[l(360/DPSI,i„)^^^j.j „ * N!UI>«,J * 3. if NFimCH - I 

lower Integer 


(Thie ie punched output.) 


It must be stated that Is the maximum possible number of 
output records which would be produced by the program version, only If 
all input options were used and only If tl.'e configuration used LLL panels 
and only if MAXALP orientations were input. 

Normally the configuration will not have LLL panels but will have 
NP panels instead. Also the number of orlencatlons used is often less 
than MAXALP and will be NALPHA instead. 

Thus to get a better estimate of the number of output records 
which will be produced when running a specific config>.s'ition, do the 
following. 

Use tho R(g( expression ^iven above, but make the following changes: 

-replace LLL value by the value of N? for the body (best estimate 
should be used here.) 

-replace MAXALP value by value of NALPHA (section 3) 

-replace MAXINF value by value of NINFLO (section 3) 

-in the R2 formula, replace the following: 

replace [NSECMA * (NIPMAX + 1)] by a value e ;ual to the 
number of cards in set 3 of the input card deck (refer 
to section 3). 

-replace NR.\D^jj by value of NRAD 
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-replace DPSljjj^jj by value of DPSI 

As an example of Che preceding descriptions, use the program listed in 
Appendix B. This program has dimensions: 

LLL - 2596, MAXALP » 6, MAXINF • 500, NSECMA - 60, NIPMAX • 45, 

NRAD - 51, and DPSI - 1.0 

r.idx din 

Using these values in the formula gives the maximum possible 

ntimber of output records by this program version: 

R^iot “ 363,782 records 

This far exceeds the 30,000 record limit imposed at Che Penn State 
computer facility. Hx^ever, a typical configuration run with this version 
of the program might consist of: 

NP * 968 panels 
NALPHA - 6 
NINFLO - 88 
NRAD » 21 
DPSI *15.0 
ITMAX - 20 
NSYMET - 0 
NCALC » 0 
fn^ING - 1 
NPOINT = 0 
NPUNCH = 1 

with 877 cards in the set 3 of the input card deck. 

Using the formula for Rtot with the above numerical substitutions, 
gives the following; 

R^gj. = 25,266 output records actually produced. 

This specific configuration could be run at Penn State since it meets 
the 30,000 record limit requirement. 

Although the program dimensions are such that it fits in the 
computer core storage and executes properly, the output records may 
be excessive if the configuration uses most of the available program 
storage. As a result, it may be necessary to make tradeoffs on the 
input to the program. These tradeoffs might involve reducing the 
number of orientations submitted at one time, or attempting to reduce 
the number of body panels. 


i, 




35 


2.11 Core Storage Requlretnents Based on Program Dimension Size 

Core storage requirements for the source statements is constant. 
However, the core storage needed to store arrays and COMMON blocks is 
a function of the dimension size variable values for the program version 
being used. The program is run in single precision. For the IBM 370 
computer, this means each number in the arrays and COl'lMON statements 
occupies 4 bytes of storage. 


The maximum storage available to users at Penn State Is 560 K bytes 
core storage (K equals 1024). However, not all of the 560 K bytes are 
available directly for program storage. Son£ of the storage is taken 
by the computer for compiler operation and is not available for program 
array storage. 

An estimate of core storage required by the computer program on 
the Penn State computer facility is obtained by the following formula 
(refer to definitions in section 2. 3) : 


S. “ (S , , + S ) 

tot obj arrays 


K bytes, where K » 1024 


where is total core storage needed by the program 

Sobj is core storage needed to store the program source 
statements 

S is core storage needed to store arrays and other variables 

in COMMON block.c,. ' 


SqJjj is a constant for this program: 

118656 „ , ^ 

= • • £ 0 24 ^ 


S is a function of the program dimension size variables: 

arrays 

^arrays “ 1 I 24 * * LLL + (NIPMAX * NSECMA * 3) 

+ (17 * MAXALP) + (6 * NIPMAX) 

+ (2 * NSECMA) + (2 * MAXINF) + 281] Kbytes 


Estimation of S^ot gives an approximate idea of how much core 
storage to request. However, the total storage requested must also 
provide enough space to satisfy S^q^ as well as storage required by 
the computer compiler. 


'intaf iir 
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For exatople, the program version listed In Appendix B is dimen- 
sioned as follows: 

LLL « 2596, MAXALP « 6, N!JECMA • 60, NIPMAX -45, and MAXINF » 500 

From the formula above compute: 

= 316.69 K bytes 

It has been found that Stot plus the compiler storage requirement 
does not exceed 560 K bytes. Thus the program version of Appendix B 
can be run in the 560 K storage region (only as a category 5 run) on 
the Penn State computer facility. 

For a second example, the same program but dimensioned to a 
smaller size (listing not included in this maniial) has been written 
using: 

LLL = 2450, MAXALP » 1, NSECMV » 50, NIPM4X = 51, and MAXINF « 500 

From the storage formula compute: 

S. . = 257.69 K bvtes 
tot 

It has been found that with this smaller size version of the 
program, the source deck must be compiled into an object deck using the 
Fortran H compiler. Optimization level 2 feature on the Penn State 
system and needs 560 K for this. However, once this object deck has 
been made, it can be executed using onJ.y 280 K storage. The 280 K 
storage required means this version of the program can be run in higher 
priority job categories at Penn State and will have faster turn-around 
time. Refer to sections 2.13 and 2.14 on compiling and using an object 
deck for the program. 

2*12 Execution Time and Object Program Considerations 

Tlie iterative solutions of hundreds of simultaneous equations by 
the computer program is the process requiring most of the total 
execution time. Tlie required execution time incre».*des rapidly as the 
number of body panels (equations) increases. System time used in 
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reading equations from the auxiliary file storage device increases as 
the number of equations increases and as the number of solution 
iterations Increases. Also more time is needed as the number of body 
orientation cases increases* and the number of survey points in the 
propeller plane Increases. 

No quantitative relationships have been made between run time and 
the above program factors. ' However, Table 2.2 lists several cases run 
on the Penn State IBM 370 computer, and these give some insight into 
the run time requirements as a function of number of equations, number 
of iterations, number of orientations, and number of survey points. 

The time limi^ imposed on runs made on the Penn State facility is 
2000 seconds. This is the maximum amount of time ever allowed for any 
run. The last run on Table 2.2 required more than 2000 seconds and 
failed to give the complete solution. 

Some obvious points can be made regarding run time. When possible, 
a configuration should be submitted using the sjTnmetric input option 
(see sections 1 and 3). This will halve the number of equations which 
must be solved simultaneously and will greatly reduce run time. It is 
likely that a configuration, run on the Appendix B size program, xjhich 
uses most available program storage (i.e. has nearly LLL equations, 
MAXALP orientations, and maximum number of propeller plane survey 
points) will require excessive run time on the Penn State computer and 
will not be completed. However, no configuration this large has been 
run and the actual time requirement in this situation is not certain. 

To keep run time within limits, it may be necessary to compromise on 
runs with configurations having very many panels. For example, reducing 
the number of orientations submitted and the number of propeller plane 
survey points may be necessary Xi?hen the configuration has many panels. 

From experience, due to the lengthy execution times required on 
typical configurations having 1000 or more panels, it is highly 
recommended that the source program deck (i.e. as ic is listed in 
Appendix B) be compiled using the Fortran H compiler, Optimization 
level » 2 feature available on the Penn State computer system. The 
resulting object program deck may be punched on cards or placed on 


Table 2.2 Comparison of computing time requirements of several run cases. 
(IBM 370 computer and IBM 3330 disk pack used.) 
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magnetic storage (BAT) files. Runs are then made by submitting and 
executing this highly efficient object deck. The resulting decrease 
in run time with the object deck has been on the order of 50 percent. 

In some configurations, use of the object deck has enabled solutions 
to be obtained whereas the same configuration run using the source 
deck would have failed to finish within the 2000 second time limit. 

Refer to sections 2.13 and 2.14 for the details of compiling the object 
deck from the source deck and submitting runs with the object deck. 

2.13 Object Program — Considerations and Procedures for Compiling It 

The procedures and discussion of this section pertain specifically 
to the Penn State University computer facility which uses an IBM 370 
computer. 

It is highly recommended that the FORTRAN source program (as it 
appears in Appendix B) be first compiled using the FORTRAN H compiler, 
under Optimization level 2. The output is a very efficient machine 
language object program deck. This object deck output should be punched 
on cards or written to the magnetic disk files (BAT files on the Penn 
State facility) . 

When submitting data and running the program, the object program 
deck is submitted directly with data. Time is not lost in compilation 
since the object program has been compiled before submission. Most 
importantly, the H, opt = 2 object program is efficient and executes 
much faster than occurs when the original source deck is submitted for 
compilation and execution. As pointed out in section 2.12, experience 
has shown that as much as 50 percent execution time savings is attained 
by using the object deck directly with the data. Time saved in using 
the object deck may allow large jobs to complete within time limits 
whereas the same case performed by submitting the source program for 
compilation and execution would have failed to complete within the 
time limits. 

The follotrfing two subsections give specific instructions for 
compiling an object deck and storing it on cards and disk files, 
respectively. 
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2.13.1 Inatructlona for Compiling thfe Object Deck and Punching It 
on Garda 


To compile the aource program using FORTRAN H, Optimization 
level 2, and punching the object deck on cards, do the following. 

Submit the following cards exactly as shown (all cards below start In 
column 1 unless otherwise noted) : 

// job card (use category 5, specify T*2000, R» 30000, S“280 K) 

// EXEC FHC, FARM. SOURCE-' DECK, OPT-2' 

//SOURCE. INPUT DD * 


. FORTRAN source deck version of the program 
. (cards as they appear in Appendix B, or the BAT 
. files containing these card images. This is done by 
/* placing cards here of the form: 

/*INCLUDE Userid. $ source file name 1 
/^INCLUDE Userid. $ source file name 2 
etc. , until all source files ave given) 

Note: Userid is the identifying number a user must have for 

using the Penn State BAT file system. 


2.13.2 Instructions for Compiling the Object Deck and Writing it 
on BAT Files 

To compile the source program using FORTRAN H, Optimization 
level » 2, and writing the object deck on BAT files, do the following. 
Submit the following cards exactly as shown (all cards below start in 
column 1 unless otherwise noted) : 

// job card (use category 5, specify T=2000, R-30000, S=560 K) 

/■'USERID card for using the Penn State BAT file system 
// EXEC FXC, FARM. SOURCE- 'DECK, OPT-2' 

//SOURCE. SJSPUNCK DD UNIT-BAT, FILES-(Sfilename 1, . . . ,$filename 4) 
//SOURCE. INPUT DD * 


FORTRAN source deck version of the program 
(cards as they appear in Appendix B, or the BAT files con- 
taining these card images. This is done by placing cards 
here of the form: /*INCLUDE Userid. Ssource file name 1 

/*INCLUDE Userid. Ssource file name 2 
etc. until all source files are given) 


/* 

Note: Sfilename 1, Sfilename 2, Sfilename 3, and Sfilename 4 are the 

names of the four BAT files which sequentially store the output 
object program deck. These names must fit the naming con>rentions 
for BAT files. (See references 6 and 8.) 
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The potential flow program requires four BAT files to store its 
object deck. Four files are required regardless of the dimension size 
of the program being compiled. 

Instructions for running the object program with data are given 
In section 2.14 of this manual. 

2 . 14 Operating Instructions — Required Job Control Cards 

The following instructions and job control cards apply specifically 
to the source program version listed in Appendix B or Its corresponding 
object deck. ITie instructions will apply to any size version of the 
program if the job control cards establishing file 9 and file TjO storage 
are modified according to the rules of section 2.9. Also the job control 
cards apply specifically to the Penn State IBM 370 computer facility. 

The program in Appendix B has been dimensioned as follows: 

(See sections 2.3 and 2.4 for definitions.) 


NSECMA 

NIPMAX 

MAXALP 

MAXINF 

LLL 

MD 

LLLHAF 
DPS I 


nan 


NKAD, 


max 


60 

45 

6 

500 

2596 

2602 

1298 

1.00 

51 


Subsections 2.14.1 and 2 .14 .2 below give operating instructions for running 
the program from the sov»‘ce deck and object deck, respectively. These 
are applicable to the computer system at Penn State. 


2.14.1 Running Program by S ubmitting Source Deck and Data 

This is the method of using the program which is simpler 
but less efficient. The source program r'uns slower, and it is not the 
recommended method. However, this may be the only way the program can 
be used at computer facilities different from that at Penn State. 
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The Job control cards given in the following allocate file 9 and 
file 50 disk storage corresponding to the above program dimension sizes. 
For using the program with different dimensions, these Job control 
cards must be mo(*lfied using the Instructions In section 2.9. 

Supply the following cards, In the order shown, to compile and 
execute the source program and process data. Punch these cards 
exactly as shown starting In card column 1 unless otherwise noted: 

// job card (use category 5, specify T'=2000, R* 30000, S*560 <) 

// EXEC FHCG 
//SYSIN DD * 


. source program: (cards as they appear in Appendix B or 

. the BAT files containing these card 

. images. This is done by placing cards 

. here of the form: 

. /^INCLUDE Userid.$source file name 1 

. /^INCLUDE Userid. $source file name 2 

. etc. until all source files are given.) 

//DATA.FT09F001 DP UNIT-SYSDA, SPACE«(CYL, (123, 1) , RISE) , 

/ / DCB- (RECFM- VS , LRECL-3116 , BLKS I2E-3120 ) 

//r>ATA.FT50F001 DD UNIT-SYSDA, SPACE-(CYL, (1,1) ,RLSE) , 

/ / DCB- (RECFM- VS , LRECL-40 , BLKS I2E-44) 

//DATA. INPUT DD * 


. input data cards punched according to format of section 3 
. of this manual. Also may have data on BAT files included 

. here by punching the following cards here of the form: 

/*INCLUDE Userid. Sdata file nama 
. Include enough of these cards to supply all necessary 
. data B.\T files. Cards and BAT file include cards may be 
. intermixed. 

• 

— This concludes the operating procedure. — 


2.14.2 Running Program by Submitting Object Deck and Data 

This second operating method is recommended. It Is more 
complex, however, because the object program deck must be obtained using 
the methods of section 2.13.. With the object deck available on cards or 
BAT files, jobs may be run using the procedure given in this subsection. 
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As In section 2.14.1, the job control cards given here establish 
file 9 and file SO storage space to aeet the needs of the program in 
Appendix B whose dimensions were listed above. These job control 
cards may be modified to suit a program of different dimensions. This 
is done by following the procedures in section 2.9. 

Supply the following cards, in the order shown, to execute the 
object program deck and process data. Punch these cards exactly as 
shown starting in card column 1 unless otherwise noted: 

// job card (use category 5, specify T «2000, R»30000, 3 *560 K) 
// EXEC FHG 
//DATA. DECK DD * 


object program: (card deck or images of cards from BAT 

files as produced by section 2,13 instruc- 
tions. To submit program from BAT files, 
place the following four cards here of the 
form: /*INCLUDE Userid. $object file name 1 

/^INCLUDE Userid. $object file name 
/^INCLUDE Userid. Bobject file name 
/^INCLUDE Userid. $object file name 
where the above file names are the same 
as those used on the cards when following 
procedures of section 2.13.2.) 


//DATA.FT09F001 DD UNIT-SYSDA,SPACE«(CYL, (123,1) ,RLSE) , 
// DC,3=(RECFM*VS,LR£CL*3116 ,BLKSIZE-3120) 
//DATA.FT50F001 DD UNIT*SYSDA, SPACE- (CYL, (1,1) , RISE) , 

/ / DCB- ( RECFM-VS , LRECL-40 , BLKS IZE-44 ) 

//DATA. INPUT DD * 


input cata cards punched according to format of section 3 
of this manual. Also may have data on BAT files included 
here by punching the following cards here of the form; 

/^INCLUDE Userid. $data file name 
Include enough of these cards to supply all necessary data 
BAT files. Cards and BAT file include cards may be 
intermixed. 


— This concludes the operating procedure. -- 

This concludes the detailed descriptions of the program. In the 
next sections of the manual are descriptions of input; data and output 
data organization. 


C'J ro 




SECTION 3 INPUT DATA DESCRIPTION 


The Input to the prog rem conalets of eight card sets, each set 
containing one or toore cards. No card set is ever omitted. Some input 
options are available and are explained below and in the card 
descriptions: 

1. Symmetric body input option: This may be used if body is 

symmetric about the y ■ 0, x-z plane. No closed separate 
bodies such as tip tanks or nacelles can be present. The 
wing, if present, must be symmetric, and no body sideslip 
angles may be specified. Then only the left side of the body 
geometry need be provided by the user in card set 3. 

2. Regular body input option: If any condition in (1) above is 

not met then the regular input option tnust be used. The user 

must provide the complete body surface geometry in card set 3. 

3. Test run option: This is used only to genera*‘c and print the 

body panel geometry. No flow calculations are made. Two 
purposes of this option are: 

a. To check for errors in the geometry input. 

b. If user is not certain which index number will be assigned 
to a certain panel by the program, this option allows the 
user to identify each panel by index nunber (as is required 
when specifying Inlet/outlet panels in card set 5). In 
using this option, all card sets are input as instructed 
except no Inlct/outlet panels may be specified (NINFL0“0 
in card set 3). After making the test run, the proper 
index numbers of inlet panels are identified from the output 
panel geometry and can be input in curd set 5. Then the 
card deck may be run as a normal flow prediction run. 

4. Printing and punched output options: Vser may select to print- 

out the panel geometry or not to print it. Note for a test run 
this geometry' is automatically printed (see output description 
in this man<ial) . Also, user may elect to punch or not punch 
certain propeller plane flow output on cards. 




'■"W 
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Thia card Input daacrlpcion^ below, la written In general fora 
applicable to any dintenalon size version of the prograa which the user 
has available. That Is, the value Halts on th« Input variables In the 
card descriptions are expressed In terms of the following prograa dlaen- 
slon size variables: NSECMA, NIPMAX, MAXAL?, MAXINF, LLL, MD, LLLKAF, 

DPSI^^j^, and defined In detail In the aaln program 

and govern the dlaenalonlng of arrays. They are explained In detail In 
section 2 of this manual on changing prograa dimensions. The definitions 
of the above size variables are: 


NSECMA - Maximum number of body cross sections which may be input. 

NIPMAX ■* Maximum nunber of periphery points allowed to be input 
around any complete cross section. 

MAXALP - Maximum number of input body orientations allowed. 

MAXINF - Maximum nurber of inlet/outlet panels which may be 
present over the entire body. 

LLL - ■(NSECMA-D’^CNIT’MAX-l), Maximum possible uunber of body 
panels on entire body. 

>fl) - -(LLL + .MAXALP) 

LLLHAF - -(LLL/2) 


DPSIg^u - Minimum allowed azimuth angle spacing for points on 
propeller plane. 

NRAO - Maximum allowed number of points along each azimuth in 
the propeller plane. 


Keep in mind Chat if the input is being used in the program version 
listed in Appendix B then the dimension variables have the following 


values : 

NSECMA • 60 

NIP!IAX • 45 

ma:cup » 6 

MAXINF » 500 

LLL ■ 2596 
MD - 2602 
LLLHAF • 1298 
DPSI„in - 1.0 

NRADmax • 51 
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Ua« any contlatcnt family of units for tha vtlocitios and coordl- 
natas throughout tha data. Dlnenslonlass quantitias are so nonad. All 
angles are in degrees. 

Integers are right justified in card cwlunn fields and are identified 
as integers in the description. 

Floating point numbers are punched anywhere in the given column 
field and must include the decimal point. These type numbers are 
identified as floating point in the card description. 


Card Set 

1: Data or Run 

Identificat 

ion (2 cards) 

Card 1.1 

Title Card 1 



Notes 

Variable 

Columns 

Description 

(1) 

SYliBOL(I) , 
1=1,80 

1-80 

First part of the Input title. 

Card 1.2 

Title Card 2 



Notes 

Variable 

Columns 

Description 

(2) 

SYMBOL(I), 

1-31,160 

1 - 8C 

Second part of the input title. 

Note;! for 

Card Set 1 




(1) Contains first part of a title containing any desired 

identifying Information in letters, numbers, or symbols. Card 
must always be prc>sent in deck even if totally blank. 


(2) Second card for continoation of the identifying information. 
Contains any letters, numbers, or symbols. Card must always 
be present in deck even if totally blank. 


Card Set 2; Free Stream Velocity — Body t')rlentatlons 

Set specifies the body orientations for which flow solutions are to be 
obtained. (Set contains 3 to MULA.L? + 2 cards.) 

Card 2.1 Free Stream Velocity 

Variable Column s Description 

V 1 - iO Floating Point. Magnitude of free 

stream velocity. (1.0 is convenient) 


Notes 

( 1 ) 
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Card 2.2 

Number of 

Orientations 

Notc0 

none 

Variable 

NALPHA 

Column 

1 

Cards 2.3 

Body Orientations (NALPHA 

Notes 

(3) 

Variable 

ALPHA(I) 

Columns 

1-10 

(4, 5) 

BETA(I) 

20-29 


D<acrlt>cion 

Integer. Nuo^er «f body orienta~ 
tions (pairs of angle of attack and 
a idee lip) this run. Value from 1 
to MAXALP. 

cards. See note (2).) 

D.fscrlptlon 

Floating point. Body an;;lc of 
attack in degrees for orientation 
case I. (Value of 0 to 360., or 
0 to -360.) 

Floating point. Body sideslip 
angle in degrees for orientation 
case I. (Value of 0 to 360., or 
0 to -360.) 


Notes for Card Set 2 

(1) All output velocity quantities are nonditocnslonalizcd with 
respect to V so the actual magnitude of V is not Important 
and a value of 1.0 is rcconsnended for convenience. 

(2) Card 2.3 is repeated NALPRiV times, each with an angle of 
attack and associated side slip angle. The first card 2.3 
is orientation case 1, the second card is case 2, etc. 

(3) ALPHA is the angle of attack between the free stream velocity 
(projected onto the reference body x-z plane) and the reference 
body X axis. It is positive for body nose upward. See 
Figure 3.1. 

(4) BETA is the angle between the free stream velocity (projected 
onto the reference body x-y plane) and the reference body 

X axis. It is positive for body nose to the right. See 
Figure 3.1. 

(5) All 3ETA(I), 1*1, M/LPHA must be zero if the s-ymmetric body 
input option (NSYME* ■ 0) is specified in card set 3 below. 




jMjiiaii 
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Card Set 3; Symmetry Option and Body Paneling Input j 

This set of cards specifies the use of symmetry or non-s>Tnmetr>' in panel | 

input. The body cross section geometry points are given here. The left I 

side of the geometry is given if symmetry is used. Otherwise both sides i 

of the geometry must be given. All geometry is in terms of the reference | 

body aircraft axis system (x foreward, y right, and z down). The origin j 

may be anywhere desired except if jyt Vvjtry is used. If symmetry is used, 
the origin must be on the plane of symmetry. See also section 1 in the 
manual which explains body paneling. (Set contains 11 to 
((NSECMA * NIPMAX) + NSECMA + 2) cards.) 

Card 3,1 Symmetric Body Input Option Choice 

Notes Variable Column Description 

(1,2,3) NSYMET 1 lnteger*0 if symmetric option used. 

Integer**! if option is not used. 

Card 3.2 Number of Body Cross Sections 

Notes Variable Columns Description 

(4) NSECTO 1-2 Integer. Total number of body 

cross section descriptions given 
below. (3 NSECTO i NSECMA) 

Card Groups 3.3 Cross Section Descriptions 

NSECTO card groups 3.3 are given in sequence. Each group is associated 
with a body cross section description. Each group is composed of one 
card 3, 3. A followed by several cards 3.3.B as shown below. Groups are 
in sequence from 1 (first and front section) to NSECTO (last section). 

Refer also to section 1 of this manual. 


Card 3. 3. A Cross Section Identifier 


Notes 

Variable 

Columns 

Description 

(5) 

NSEC 

1-2 

Integer. Cross section .sequence 
number. (1 < NSEC < NSECTO) 

(6) 

NIP or 
NCOUNT(I) 

10 - 11 

Integer. Number of periphery 
points on this I th = NSEC th 
cross section. 

(3 <NIP < [NIFMLX+11/2) if 
NSYMET=0 

(4 <NIP ::^NIPJi2lX) if NSYMET==1 

(7) 

NEND or 
NFLAG(I) 

30 

Integer. Section repeat signal 
for this I th = NSEC th section. 
Equals 0 or 1 as explained in 


note. 


50 


Cards 3>3>B Section Periphery Points (NIF of these cards given for 
section NSEC identified in card 3. 3. A above.) 


Notes 

Variable 

Columns 

(8) 

X or 
PX(I.J) 

1-10 

(8) 

Y of 
PY(I.J) 

20-29 

(8) 

Z or 
PZ(I,J) 

40-49 


Description 

Floating point coordinates of the 
J th point around the periphery 
of the 1 th « NSEC th cross sec‘- 
tion. Coordinates relative to 
the reference aircraft coordinate 
system. 


Notes for Card Set 3 


(1) Symmetric input option (NSY>IET«0) should be used if possible 
to reduce solution time and cost. 

(2) Symmetry option (NSY>tET=0) is possible only if all the 
following rules are met. 

- All BETA (cards 2.3) must be zero. 

- Origin of the body reference aircraft coordinate system 
must be on the body plane of s>Tnmetry with the plane 
containing the x and z axes. 

- Configuration may not have discrete closed bodies to the 
left or right of the plane of symmetry even if such 
bodies have a counterpart on the other side of the plane 
of symmetry. Example is engine nacelles on the wings. So 
twin engine aircraft cannot use the symmetry option. 

- Only the inlet/outlet panels on the left (-y) side are 
specified in card set 5 and such panels have symmetric 
counterparts on the right side of the body which will be 
generated automatically by the program. 

- If a wing is to be modeled (card set 6 below) the wing must 
be symmetric with the wing root located on the body plane 
of symmetry. 

If any rule above is not obeyed then must specify NSYMET=1. 

(3) If the symmetry criteria are met and NSYMET-0 is used then 
read in only the left (-y) side of the body geometry points 
on cards 3.3.B. See Figure 3. 2. a. 

If the general non-symmetric input option (NSYMET=1) is used 
then input all points on a cross section on cards 3.3.B. 

See Figure 3.2.b. 

(4) A body cross section is at a constant x or nearly constant x 
location on the body, A physical cross section may be 
repeated with a new point description, but the repeat Is 
counted as an additional section and contributes to the 
NSECTO amount of sections. 



input option 


J • : 
l<3 1 2 
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(5) Crosi? sections are input sequentially so NSEC must be con- 
secutively Increasing each time a card 3. 3. A is punched. 

The firot NSEC must 1 then 2, 3 ... NSECTO. 

(6) NIP value is read from cards but renamed and inserted In 
array entry NCOUNT(NSEC) . If symmetry option is not used 
(NSYMET “ 1) then NIP is the m^er of points completely 
around the section and the first (top) point Is repeated as 
the last NIP th point. See figure 3.2.b. Refer also to 
section 1.2.1. 

If symmetry input option is used (NSYMET ■ 0) , then NIP is 
the numbe'* of points around the left half of the cross section 
with the first point on the top centerline and the NIP th 
point on the lower centerline. See Figure 3.2. a. 

NIP on section NSEC must equal NIP on adjacent sections 
between which panels will be generated. 

(7) NEND is read in but program stores it under a new name in 
array entry NFLAG(NSEC) . Use NEND « 0 if this section NSEC 

is not going to be repeated with a second description numbered 
NSEC + 1. 

Or use NEND * 0 if this is the last input section of the 
entire configuration. (NEND = 0 signals the program to 
generate panels between sections NSEC and NSEC +1.) See 
section 1.2.1 of this manual. 

Use NEND * 1 if this section NSEC, having NIP points, will 
be repeated as section NSEC + 1 with a new description using 
a different number of periphery points or the same nuniber of, 
but different, points. 

Or use NEND = 1 if the following applies; This section NSEC 
is the last input section of a discrete body in a multi- 
body configuration and another body follows this one. The 
next section NSEC + 1 is the start of a new discrete body panel 
network. For example, NSEC is the last section on the left 
nacelle (NFL.\G = 1 for this) , and section NSEC + 1 is the 
first section at fuselage nose. (NEND = 1 signals program 
not to generate panels between sections NSEC and NSEC +1.) 

NEND = 1 can never appear in two consecutive cross section 
inputs . 

(8) The coordinates X, Y, and Z are read from input and the pro- 
gram saves them under new names re 3 p.i!ctively in the array 
entries PX(NSEC,J), PY(NSEC,J), and PZ(NSEC,J) where J is the 
J th point around section NSEC as by the following rules. 






The cards 3.3.3 (points X,Y,2) must be input in the sequential 
order illustrated by Figures 3 *2. a and 3*2.b: 

- For NSYMET * 0 the first card must be the top centerline 
point. Points are input consecutively clockwise around 
the left side of the section. The last card of the section 
is the NIP th point and it must be on the lower centerline 
of the section. 

- For NSYMET « 1 the first card must be the upper point on 
the section. Points are input consecutively clockwise 
around the left side to the right side of the section. The 
last card is the NIP th point and must be an exact repeat 
of the first po*^nt. 

Card Set 4: Test Run Option and Panel Geometry Print Options (Set 

always contains 2 cards.) 

Card 4.1 Panel Geometry Print Option 


Notes 

Variable 

Column 

Description 

(1) 

NLIST 

1 

Integer » 0 if printout of panel 
geometry table is desired. 
Integer = 1 if printout is not 
desired. 

Card 4.2 

Geometry Check 

Run Option 


Notes 

Variable 

Column 

Description 

(2) 

NCALC 

1 

Integer = 0 if this is normal run 
for complete flow solutions. 
Integer = 1 if this is a test run 
to check generated body panel 
geometry. 

Notes for 

Card Set 4 




(1) Generated body panel geometr>" table printed UvSing NLIST ~ 0 
contains the following: 

^ Body panel index numbers (in sequence generated). 

-- Input corner points used to make panel. 

“ Computed control point coordinates of all panels. 

- Surface area of each panel. 

(2) Check run option (NCALC = 1) is used to get a listing of 
generated body panels. The listing helps check user 
geometry input and enables user to identify the index 
numbers of panels he wishes to specify as inlet or outlet 
panels if these are not known. 

NCALC = 1 stops program after panel generation and does not 
allow flow solution to continue. 
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NCAIC <■ 1 overrides the NLIST option and always causes the 
geometry to be printed. 


Use of NCALC " 1 assumes that user does not know the Index 
numbers of any panels which must be specified as Inlet or 
outlet panels for the flow solution. Therefore the data deck 
must be organized as follows If a test run Is to be made. 


- Set NCALC ■ 1 on card 4.2. 

- Specify NlNiT-O « 0 on card 5.1 below. 

- All other card sets contain data as if this were a normal 
run. 


After making the test run, the geometry output will allow the 
Inlet and outlet panels to be identified by index number. 


- Modify Card Set 5 by specifying the Inlet and outlet panels 
as found by the test run. 

- Change to NCALC » 0 on card 4.2. 

- Resubmit data as a normal run to get flow solutions. 


Card Set 5: Inlet and Outlet Panel Specifications (Set contains 1 to 


MAXINF + 1 cards.) 


Card 5.1 Number of Inlet and Outlet Panels 


Notes 


( 1 ) 


Variable 


NINFLO 


Columns 


1-4 


Description 


Integer. Number of inlet and out 
let body panels. Specified on 
cards 5.2 below. 

(0 i NINFLO i MAXINF) If NSYMET «= 1 
ro ^ NINFLO < (MAXINF/2) ] if 
NSYMET = 0 


Cards 5.2 Inlet/Outlet Panel Index Numbers — Velocity Ratios (Omit this 
card of NINFLO = 0, othewise provide NINFLO of these cards.) 


Notes 

Variable 

Columns 

Description 

(2) 

INDEX(I) 

1-4 

Integer. Panel index number of 
the I th inlet or outlet panel as 
generated by program. 

(3) 

FLRAIO(I) 

20-29 

Floating point. Inlet or outlet 
velocity ratio on panel INDEX(I). 


= 0.0 if panel is solid wall (not 
an inlet or outlet) . 

> 0.0 if inflow velocity (INDEX(I) 
is an inlet panel). 

< 0.0 if outflow velocity (INDEX(I) 
is an outlet panel) . 




n^i.r I 
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Notes for Card Set 5 

(1) For syniaetric option (NSYMET ■ 0) , NINFLO is the number of 
inlet and outlet panels only on the left (-y) side of the 
body. 

If the symmetry option is not used (NSYMET ■ 1) , then NINFLO 
is the total number of inlet and outlet panels on the entire 
configuration. 

(2) If NSYMET - 0, then the INDEX(I) of the left (-y) side 
Inlet/outlet panels are specified. The ones on the right 
side of the body are automatically taken care of by the 
program. 

If NSYMET ■ 1, then the INDEX(I) values of all inlet/outlet 
panels of the configuration must be given by the user. 

In either of the above situations, the INDEX(I) values on 

cards 3.2 may be input in any order desired and do not have 

to be in any special sequence. 

If user does not know which panel index nuniers correspond 
to panels which are on inlets and outlets (as may occur with 
a new untried input geometry in card set 3), it will be 
necessary to make a test run to get a listing of generated 
pane's. To do this set NCALC ■ 1 on card 4.2, set NINFLO to 
0, and follow the test run procedures given in note (2) at 
the end of card set 4. 

(3) FLSATO is a non-dimensional velocity (velocity/V) which is 

specified normal to the panel surface along the inward normal 

(positive value) or along the outward normal (negative value). 

The solid wall boundary condition will be relaxed in the 
equation for this panel nua^er IND£X(I). 

If NSYMET - 0, each FLRATO(I) value will correspond to left 
(-y) side panels INDEX(I). Program automatically assigns 
the same FLRATO(I) values to the right side panels. 

FLRATO(I) « 0.0 is valid but should not be used. It is better 
to delete this panel INDEX(I) from the card set 5. 


Card Set 6; Wing Input Data (Sec contains 1 to MAXALP + 2 cards.) 

Card 6,1 Wing Input Option 

Notes Variable Column Description 

none NWING 1 Integer » 0 if wing is not modeled. 

Integer • 1 if wing is modeled as 
given in cards 6.2 and 6.3 below. 
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Card 6.2 

Wing Geometry 

(Delete ( 

this card if NWING - 0) 

Notes 

Variable 

Columns 

Description 

(1) 

XQR 

1-10 

Floating point coordinates of the 

(1) 

YQR 

11 - 20 

wing root quarter chord location. 

(1) 

ZQR 

21 - 30 

In the reference aircraft coordi- 
nate system. 

(2) 

SPAN 

31 - 40 

Floating point. Wing span as seen 
in planform view. 

(2) 

CHORD 

41 - 50 

Floating point. Wing root chord 
length. 

(2) 

DIHED 

51 - 60 

Floating point. Wing dihedral 
angle in degrees. Positive for 
wing tips upward from horizontal. 

(2) 

SWEEP 

61 - 70 

Floating point. Sweep angle of 
wing quarter chord line in degrees. 
Positive for aft sweep. 

Card 6 . 3 

Wing Lift Coef 

ficients 

(Delete these cards if NWING » 0) 

(If NIXING 

* 1, repeat card 6.3 NALPHA times.) 

Notes 

Variable 

Columns 

Description 

(3) 

Notes for 

CL(I) 

Card Set 6 

1-10 

Floating point. Value of wing 
lift coefficient corresponding to 
ALPHA(I) of the I th orientation 
case on card 2.3. 


(1) Wing root is the location midway between wing tips; it is 
not the location of the wing - fuselage interface. See 
figure 3.3. 

If the symmetric body input option (NSYMET ■ 0) has been used 
in card set 3, then the following must be met for the wing: 

-- The wing must be symmetric about the body plane of symmetry, 
(i.e. YQR » 0.0 is required.) 

(2) See figure 3.3 for definition of angles. 

(3) Put one CL value on each of the NALPiL\ cards 6.3. These cards 
must be in sequence corresponding to cards 2.3 such that the 
first card 6.3 has CL(1) corresponding to ALPHA(l) on the first 
card 2.3 ... etc. 

CL(I) » 0.0 is valid and has the same effect as not modeling 
the wing for orientation case I. 
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Card Set 7; Propeller Plane Description (Set has either 1 or 4 cards.) 


Card 7.1 Propeller Plane Specification Option 


Notes 

none 


Variable 

NPOINT 


Colucm 

1 


Description 


none NPOINT 1 Integer ■■ 0 If no propeller plane 

is defined. 

Integer ^ 1 if propeller plane is 
modeled as defined on cards 7.2 and 
7.3. Propeller plane flowfield 
predictions will be made. 

Card 7.2 Propeller Plane Location and Orientation (Delete this card if 


Description 


Floating point coordinates of the 
hub or center point location of 
the propeller plane. In the 
reference aircraft coordinate 
system. 

Floating point. Built-in vertical 
tilt angle of propeller plane in 
degrees. (Value of 0 to 360. or 
0 to -360.) 

Floating point. Built-in sidward 
tilt angle of propeller plane in 
degrees, (Value of 0 to 360. or 
0 to -360.) 

Floating point. Reference propel- 
ler plane radius. (Propeller 
blade length.) 


NPOINT 

Notes 

- 0) 

Variable 

Columns 

(1) 

XHUB 

1-10 

(1) 

YHUB 

11 - 20 

(1) 

ZKUB 

21 - 30 

(2) 

.UPHP 

31 - 40 


BETAP 


RADIUS 


41 - 50 


51 - 60 


Card 7.3 Distribution of Flow Prediction Points on Propeller Plane 


(Delete this card if NPOINT - 0) 


Notes 

(5.7) 


Variable 

DPSI 


(6,7) 


DRADUS 


Columns Description 

1-10 Floating point. Azimuth angle 

increment (degrees) for spacing of 
radial rows of point. s around the 
propeller plane. Value must be 
positive in the range from DPSInin 
to 360.) 

11 - 20 floating point. Fractional per- 

centage of RADIUS (card 7.2) giving 
the radial spacing of points along 
a constant azimuth in the propeller 
plane. Dimensionless value of 0.0 
or any positive number. 
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Card 7.3 (continued) 


Notes 

Variable 

Columns 

Description 

(6,7) 

NRAD 

21 - 22 

Integer. Number of points to b< 
distributed radially along each 
azimuth position. Value from 1 

CO NRADgiax* 


Card 7.4 Propeller Pland Flow Field Punched Card Output Option (Delete 
this card if NPOINT - 0) 

Notes Variable Column Description 

(7) NPUNQI 1 Integer *0 if no punched card 

output desired. 

Integer ■ 1 if desire punched 
card output. 

Notes for Card Set 7 

(1) Propeller plane aay be positioned at any location desired. 
(See figure 3.1.) It need not be positioned on the plane of 
syometry, cvon if NSY.MET - D. (i.e., if NSYMET « 0, 

YHUB ■ 0.0 is not required.) 

Also the hub may be located within the body. For e:<ample, 
the hub inside of a spinner may be modeled. However, flow 
predicted at propeller plane regions inside the body will 
have no physical meaning. 

(2) ALPHP is a fixed angle of pitch between the propeller thrust 
axis and the body reference x axis. It is the angle as seen 
projected on the body x-2 plane, (i.e. as seen in the side 
view of the body drawing.) Angle is positive if propeller 
thrust axis is tilted upward. See figure 3.1. 

(3) BETAP is a fixed yaw angle between the propeller thrust axis 
and the body reference x axis. It is the angle as scon 
projected on the body x-y plane, (i.e., as seen in the top 
view of the body drawing.) A yawed propeller plane is per- 
mitted even when the sjannetry option has been specified, 
(i.e., if NSYlffiT - 0, BETA? need not be zero.) BET.\P is 
positive for thrust axis tilt toward the right (■»•;.') side of 
the body. See figure 3.1. 

(4) R.\DIUS is used as a scale length for nondimen.sionalisiny, the 
radial positions of points in the propeller plane. Also the 
RADIUS length determines how much propeller plane area is 
mapied in conjunction with values of DPvADUS and NRAD given 
in card 7.3. Choice of larger R<\DIUS will move propeller 
points farther out from the hub. See figure 3.4. 




PSI • 0 



(a) INRAD-1]*IDRADUS3 ■ 1 


PSI - 0 



i.b) INRAD-1]*(DPJU)US] - A, less than 1 

PSI - 0 



(c) [NBAD-lJ’OfDRADUS] • A, acre than 1 


Figure 3.4 Distribution of survey point:? on propeller 
plane viewed in thrust direction. 


i 


(5) DPSl should be chosen such thee ic divides 360.0 Into whole 
percs. (l.e. (360/DPSI) ■ Integer). This ensures an even 
azimuthal spacing or mapping of points around Che propeller 
plane. See figure 3.4*. 

(6) NRAD always includes Che hub point. The flow at the hub will 
be calculated repeatedly as each azimuthal row of points Is 
examined. 

The choices of NRAO and DRADUS, In combination, will determine 
Che radial spread of points over the propeller plane of RADIUS 
size. The point spread la determined by one of the following: 
(see also figure 3.4) 

- if t(NRAD - l)*DRADUSl - A, less than 1: Points will be 

spread evenly from hub to a distance A*^RADtUS, inside of 
the RADIUS circle (Figure 3.4.b) 

- if ((NRAD - 1)*DRADUS) • A, more than 1: Points will be 

spread evenly from hub to a distance A^RADIUS, outside of 
the RADIUS circle (Figure 3.4.c) 

- If f(NRAD - 1)*DRADUS] • 1.0: Points will be spread evenly 

from hub to the edge of the RADIUS circle. (Figure 3. 4. a) 

(7) If card punched output Is made (NPUNCH • 1), and if the 
punched output is to be input lO the propeller analysis pro- 
gram of reference 3, Chen the following conditions must all 
be met: 

- DPSI must be set equal Co 13.0, 

- DRADUS must be set equal to 0.05, and 

- NRAD must be set equal Co 21 


Card Sec 8: Solution Iteration Control Card (1 card) 


Variable 

ITJIAX 


Columns 

1-2 


( 1 , 2 ) 


li - 20 


Description 


Integer. Maximum allowed iter.i- 
tions in Gauss Seidel Solution for 
source strengch.H. Value from 1 to 
99, but recommended values are 13 
to 20. 

Floating point. E.xponent on the 
value of the largest permitted 
relative error. Defines solution 
convergence limit. Recommended 
value is -4.00. 
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Notes for Card Set 3 r 

I 

(1) Target solution error residual is given as: • 

error = 10.0 Thus ERR should be a negative number to get \ 

useful solutions. ! 

i 

(2) All body orientation cases, (ALPHA(I), BETA(I)], for 1*1 \ 

to NALPhA, are solved together during each iteration. All 
orientation cases must converge to the error limit within 

ITMAX iterations. If one or more cases fail to converge, all 
solutions which did conver-'® will also be lost. Program will 
abort and no solutions will given. 


This concludes the input data descriptions. 


Appendix C describes a sample run case and shows an example input 
data deck. 




SECTION 4 OUTPUT DATA DESCRIPTION 


All output is written to a standard 132 character-per-line printer* 
Additionally, if NPUNCH » 1 is specified, certain parts of the propeller 
plane flow field prediction output is punched on standard 80 column 
punch cards* These cards are intended for use as input to the propeller 
aerodynamic analysis computer program described in references 2 and 3» 

The total number of printed records depends primarily upon the 
number of panels, NF, on the configuration and on the nutii^er of tody 
oricMUaLions run. Refer to section 2.10 of this manual for output 
record requirements. 

The number of punched output records depends primarily upon the 
number of points specified on the propeller plane, and the number of 
body orientations run. Again, see section 2.10. 

In the remainder of this section, the printed output is described 
followed by a detailed description of the contents and format of the 
punched card output data. 


4 . 1 Printed Ouip\it Description 

Printed output can be catogorined into six major parts. Depending 
on the choice of print options and configuration, some of those parts 
may not always be printed. The conditions required for printing of 
a part are given with the description. The following six output parts 
are described in the soiiucnce in which they appear on the printout. 


Part I : Inpu t Data Card Listing (Always Printed) 


Tb. i s 
road in. 
listing is 
detectable 
data is ]ir 
point. Liie 
abort. It 


Cable lists all of the input data card variables in the order 
The values are identified by t!ieir variable name. This input 
printed durii'ig execution of subroutine INPUT. If any 
error is diseoverod in reading a data card, that incorrect 
intod witii an error message. The. output will stop at that 
remaining data will not ho evaluated, and tiie program will 
ail data is read succossf uii>% tiiis is indicated at the end 


of the table. 






64 

Part 2; Wing Horseshoe Vortex Strengths and Geometry Table (Printed 
only if NI«NG - 1) 

This table is printed by subroutine WGEOM only if a wing has been 
modeled. The physical wing input parameters are printed first. Next 
the geometry of the wing bound vortex, along the quarter chord, is 
printed. This Is constant regardless of body orientation. Lastly 
the table gives the geometry of the two trailing vortex filaments. 

This geometry varies with body orientation, so each orientation is 
given along with the input body angle of attack, wing lift coefficient, 
and confuted strength of the horseshoe vortex. 

Part 3: Listing of Generated Body Panel Geometry' (See note 1 below. ) 
(Printed only if test run, NCALC = 1 with NLIST = 0 or 1, has been made. 
Or printed if a normal run, NCALC = 0 with NLIST = 1, has been made.) 

This tabular listing gives the geometry of the body paneling 

produced by subroutine PANEL from the input cross section geometry. 

The first column gives the panel index number. Knowledge of these 

index numbers is needed for specifying inlet and outlet panels on 

input card set 5. The next columns give the four input corner points 

of the panel. Triangular panels may be recognized by a repetition of 

one of the corner points. Next is printed the computed control point 

coordinates of the panel, relative to the body reference axes. Lastly, 

is printed on the surface area of the panel. The end of this table is 

denoted by a completion message and a tally of the total number of 

panels, NP, on the configuration. 

Note 1: For a test run (NCALC = 1) part 3 above is the last out- 

put data obtained. 

The following printed parts and punched data are obtained if a 
normal run (NCilLC = 0) was made: 

Part 4; Solution of Source Strengths, SIGMA, Table (.-Uwavs printed v:hen 
NCALC = 0) 

This table gives a summary of the progress during the iterative 
solution of equations for the panel source strengths. The iteration 
and target error limits are given. Then for each iteration performed. 
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the body orientation case, body angle of attack, and sideslip angle 
are printed along with the resultant solution error for that lteratlon4 
As an orientation case converges, this Is Indicated by a message. 

If all orientations converged within the iteration limits, this 
Is stated in a message at the end of this output. Then the remaining 
parts 5 and 6 will be printed. 


If one or more orientations fail to converge, an error message Is 
printed. All converged solutions are lost and the program Is aborted. 


Part 5; Body Surface Velocities and Pressures Table (See note 2 after 
part 6 description.) (Always printed when NCALC ■ 0) 


One of these tables is produced for each body orientation case. 

This output is printed by subroutine VELOCI. The table gives the orien- 
tation case number, body angle of attack, and sideslip angle. Free 
stream velocity magnitude and components are printed as well as the 
total number of body panels. Each panel index is listed with its control 
point. The nondimens ionali zed cartesian velocity components computed at 
the control point are printed. These components are relative to the 
body reference coordinates system. Lastly the pressure coefficient, 

CP, at the panel is given. 




Part 6; Propeller Plane Flow Field Output Table (See note 2 belot^.) 
(Printed when both NCALC = 0 and NPOI.NT * 1) 

This output gives the flow field predictions in the propeller 
plane for one body orientation case. Part 6 is divided into two sub- 
parts 6. A and 6. B described below. Part 6 output is produced by sub- 
routine VPROPS. 

Sub part 6. A Propeller Plane Output (Part 1) 

This table prints the orientation case number as well as body 
angle of attack, ALPHA, and sideslip, BETA. Free stream velocity 
magnitude, V, and components are given along with the input propeller 
reference radius, RADIUS. Refer to figure 4.1. 


.i 


The table proceeds point-by-point in the plane. The point location 
in cylindrical coordinates (R, PSI) and the cartesian coordinates of the 
location, in the body reference coordinate system, are printed. Lastly, 
the nondimensionalized cartesian velocity components, VXprop, V’Yprop and 



velociclea 
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VZpj.^p, at the propeller plane point are given. These velocities are 
relative to the body reference coordinate system. They represent the 
contribution of velocities induced by solid panels. Inlet/outlet panels, 
the wing, and the free stream velocity. Again refer to figure 4.1. 

Subpart 6.B Propeller Plana Output (Part 2) 

In this table, each 'point in the propeller plane Is Identified 
by its radial'>-azimuchal position only. The nondimensionalized axial, 
tangential, and radial induced velocity components are printed. These 
cos^onents are relative to the cylindrical coordinate system attached 
to the propeller plane (see Figure 4.1). As in subpart 6. A these 
v'elocicles represent the contributions of solid panels, inlet/outlet 
panels, the wing, and the free stream velocity. Lastly, the flow 
angles corresponding to the velocity components at a point are printed. 
These angles are rotational flow angle, outflow angle, upwash angle, 
and sidewash angle. Also these angles are named ROTANG, OUTANG, 

UPVJASH, and SIDWSH, respectively, in the program and are shown on 
Figure 4.1. See also reference 1 for the definition of these angles. 

Note 2: Output parts 5 and 6 are printed in pairs for a single 

body orientation case. The first pair of tables is printed for orienta- 
tion ease I followed by a second pair for orientation case 2, etc. 

This concludes the description of the printed output. A sample 
of printed output for an example case is shox^n in Appendix C. 


4 . 2 Punched Output Description and Format 

This punched output card data is intended for use as part of the 
input cards to the propeller performance prediction program of reference 
3. Tnese cards may be directly used in the reference 3 program xcithout 
modification. 

Punched output on 80 column cards is obtained only if NPUNCH * 1 
with NCALC * 0. This data is collected on auxiliary file 50 in sub- 
routine VPROPS. Then the contents of file 50 are punched on the cards. 
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For each body orientation case a separate group of output cards 
la punched. The output catd groups are punched In the' order natchlng 
the Input order of the body orlentatlona (input cards 2.3). 1herefore« 
there 9111 be MALPllA card groups punched. The contents of a Card 
group is given la detail r below. Mote a single ^ c 

represents the cooplete card deck needed for input data set III of the 
propeller perf oroance prediction program In reference' 3* See this 
.reference -for^ details 

Columns Indicates the card column field in which a value is 
punched. F format represents floating point nunh era. I format denotes 
integer values which are always punched right justified in their column 
fields. 

Each card group, 1: 1*1 to NALPHA, contains three se^ (types) 

of cards as follows: 


Card Set 1:. Title Cards (always 2 cards) 
Notes Variable Columns Format 


Description 


none SYJfflOL(I), 1-80 Alpha- Input title supplied by user 

I «> 1, 160 numeric on cards 1.1 and 1.2 of input 

data. 


Card Set 2; Overall Propeller Plane Orientation for Group I (1 card) 


Notes 

Variable 

Columns 

Format 

Description 

(1) 

NPSIS 

1-5 

15 

Number of azimuth positions. 

(2) 

NSAD 

11 - 15 

15 

Number of radial positions. 

(3) 

RADIUS 

21-30 

F10.4 

Reference propeller plane 
radius. 

(4) 

ATOT 

31 - 40 

F10.4 

Total propeller angle of 
attack. (Degrees) 

(5) 

BTOT 

41 - 50 

F10.4 

Total propeller plane side- 
slip angle. (Degrees) 

(6) 

ALPHA(I) 

51 - 60 

F10.4 

Body angle of attack for 
orientation case I. (Degrees) 

(7) 

BETA(I) 

61 - 70 

F10.4 

Body sideslip angle for 
orientation case I. (Degrees) 



*■ 


Hot«a for Card Sot 2 

(1) Tlilo lo eolodotod by ycogros ftte 


(2) 'llli^ lo tb« oiBii volvli-oybd^llilf^Of; 




tigu«o'':4viK^' : 


7*2 <dce 


'„•, >■ :’.::rL'.^' 

■.jM-.'--. 

' ■•• «i ■, 


for output, cord group I (body orleototteo cooo I) » AtOX ii 
^ the oet pirepollOr plane angle of Attobk betwe^ the ptopellt^ 
thruat axla and the free •trean veloelty aa eeen proiaetad ; ^ 
onto the body je*a plane* it la given by the Input an^ 

ATOT « A2jPBA(t).> AX^ (iiM 

where ALPHA(X) vna given on thei i th card 2« 3 of Input* 

AtfV is the angle given (»i lin>ut card 7*2. 

(5) For output card group I (body orientatlMi case I)« BIOT la 
the net propeller plane aideallp angle between the propeller 
thruat axis and the free atreaa velocity aa seen projected 
on the body jc^y . plane. It la. given by the Input anglea: . 

BTOT - BETA(I) + BETAP (see Figures 3.1.b and 4.1) 

was given pn the X th card 2.3 of Input. 

BETA? Is the angle given on input card 7.2.^ 

(6) ALPHA (I) Is the same value specified on Input cards 2.3. 

(7) BETA(I) Is the aaine value specified on input cards 2.3. 


Card Set 3; Propeller Plane Point— Velocity Component Cards (For each 
card group there are NPSIS * NRAD of these cards pimched.) 


Notes 

Variable 

Columns 

Format 

Description 

(1) 

PSI 

1 - 10 

F10.4 

Azimuth position of point. 
(Degrees) 

(2) 

RRATIO 

21 - 30 

F10.4 

Radial position of point, 
(Nondlmensionalized} 

(3) 

VA.XIAL 

41 - 50 

F10.4 

Axial component of velocity 
at point. (Nondlmensionalized) 

(4) 

VTANG 

61 - 70 

F10.4 

Tangential component of 


velocity at point. 
(Mondlmens ionall zed ) 


On« card for Mch point Is punchod* The card order Is such that 8BATI0 
increases aost rapidly: 

For exanple**«Aa8unis the points are spaced at every 15 degrees of 

esisuth with 11 points radially spaoed at 10 percent 

' ' 

' The card set- 3 : 0 rder will- 

first NRAD cards at PSI « 0.00 with HRATIO • p.O to 1.0 • 

second NRAD cards at PSI “15.0 with RRAITO ■ ,0.0 to 1»0 

NPSIS thNiUU) cards at RSI • 345*0 with RRATIO - 0.0 to 1.0 
Notes for Card Set 3 

(1) Axinuth angles are clockwise Increasing around the propeller 
plane disk as viewed In the direction of thrust. The first 
azimuth is always 0.0. Each azimith position is spaced DPSI 
degrees from the others. DPSI was input on data card 7.3 
(see Figures 3.4 and 4.1). 

(2) RRATIO is the dimensionless distance outward from the 
propeller plane hub to the point. It is a fraction of the 
reference propeller plane radius as follows: 

RRATIO ■ R/RADIUS 

where R is the dimensional radial distance from hub to 
point (Figure 4.1). 

(3) VAXIAL is nondimenslonalized by dividing by the free stream 
velocity magnitude^ V (on input card 2.1). VAXIAL is positive 
in the direction of propeller thrust (right hand propeller 
rotation assumed). (See Figure 4.1) 

(4) VTANG is nondimenslonalized by dividing by the free stream 
velocity magnitude, V (given on input card 2.1). VTANG is 
positive if it follows the local direction of blade rotation 
(right hand propeller rotation assumed). A positive value of 
VTANG, as seen by an observer positioned on a blade section, 
would subtract from or reduce the net onset velocity the blade 
section would see (see Figure 4.1). 

This concludes the description of the program output data. 
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AFFBHDIX A.. ■ 

Mtto Proaran md Subroutln* Ptscrlptloas 

(Milii progru glvcii first followsd by subroutiUisi la slphsbstlesl ordsr) 

Each of tbs following dsserlptlons stBamritts tha purpoast aathod, 
aad uta of tbs aubtoutlnaa. Also givao is tba nans of eosnon blocks 
rsquirad in aach. Tha principal input variablas ara llatad, ganarally^ 
in order of first appaarance. Slailarly* tha output yarlablaa ara 
llatad.'\";. 

Tha parta of the con^utar program dascrlbad ara: 

MAIN PROGRAM 
Subroutine ANGLES 
Subroutine OOEPZC 
Subroutine COFSTM 
Subroutine EULER 
Sifjroutina INPUT 
Subroutine PANEL 
Si^routlna SOLSTM 
Sv6routine SOLVE 
Subroutine VCOMP 
Subroutine VELOCI 
Subroutine VORTEX 
Subroutine VPROPS 
Subroutine WGEOM 
Subroutine WINGV 
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MAIN PROGRAM 

Purpoflft: To ostabllab progtaa array dlMnatoo alsaa and dafina 

valuaa of axacutloa tina dlMaaipning variablaa. Alao 
to link aiib rout Inaa toga thar and anaura thay ara 
^ axacutad, In tha propar aaquanea. 

Matbod: Daar apaelfiaa thi valuaa of four dlaanaion aiaa control 

variablaa and proparly dloanalona arraya in tha dioMn- 
alon atataaant at tha atart of tha program, hafar to 
aaetlon 2 of thia manual for dataila. 

•* V . ■ : - . if 

4'’- ■ - i/ 

Tha main program tikaa tha uaar auppllad airing vari- 
ablaa wd from thaae, caleulataa othar airing variablaa. 
All of thaaa quantitiaa are chaekad to anaura they hava 
valid valuaa aa dafinad by tha rulaa in aaetlon 2 of 
thia manual. 

If all tha dimanaion airing variablaa ara valid* the 
main program exeeutaa tha aubroutlnea and uaaa Che 
variablaa for automatic axacutlon time dimanaioning 
of all aubrouclna dimanaion atatamanta. 

Arraya in the COMMON atatamanta of aubroutlnea ara not 
autoeiacically dimanaionad by the main program. Thaaa 
muaC be dimanaionad by uaar by actually changing the 
COMMON carda in tha at^routinaa if a change in program 
dimanaioning variablaa ia made. Sea aaetlon 2 of thia 
manual. 

Finally* if any of the dimanaion alre variable values 
la invalid* the main program printa an error message 
identifying the problem. Program then aborts. 

Common Blocks 

Required: INPUTS, SYMTRE* WINGl 

Input : (User ch«tgaa tha following cards if program 

dimensions are to be modified). 

NSECMA Maximum limit on number of body cross section 
daaeriptlona permitted in defining geometry. 

NIPMAX Maximum limit on number of points specified 
around any complete body cross section. 

MAXALP Maximum number of body orientations which may 
be solved in one run. 

MAXINF Maximum total number of inlet and outlet panels 
which may be specified over the eurire body. 


I 


! . 


i 




] 

1 


i 


I 


I 








Subroutine* 

Called: 


Error Stops: 
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Other Input : (obtained from data input) 

NWXNG 


Wing input option ■ 0 if no wing modeled* 
* 1 if wing is modeled. 


NSYMET 


Symmetric body input option 

■ 0 if synnetry used. 

■ 1 if synnetry not used. 


NALPHA tiuober of body orientations input. 


Output ; (Dimension size variables generated by program 
based on user supplied ones) 


LLL 


liaximum total niimber of body panels which can 
possibly be generated consistent with geometry 
input limits of NSECMA and MIPMAX. 


MD 


Sum of LLL and MAXALP, will equal nue^er of 
coefficients in each equation for unknown 
source strengths. 


LLLHAF LLL/2, represents half the total maximum number 
of body panels which may be generated. 


Other Output : 


ISTOP Dimensioning variable error detect signal 

■ 0 if no such errors found. 

■ 1 if one or more errors detected. 


INPUT, VCOMP, WCEOM (if NWINC *• 1), PANEL, COEFIC (if 
NSYMET - 1), SOLVE (if NSYMET - 1). COFSYM (if NSYMET-0), 
SOL.SYM (if NSYMET - 0), VELOCI, and VPROPS 


If ISTOP • 1 

Print error message about dimensioning variables, run 
stopped. 
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SiitorotttiM AIK&tS 




Hidtedt 


; To flov angl** outflow angA«, 

upfwttti v^lOt and oldiwatb anglo «t « point on to* 
■^wpollf? plnni- 

Tho vototlodnl flow «nd mitflow onglM aro eoi^utod 
In th« propollot piano eyllndrietl eootdinata ayotan 
00 doflnod in tafwranoa 1 and Flguro 4.1 of thia annual. 
Opwaah and aldawdab angina ara oot^utad la tha rafaranea 
body coordinOta ayacaa an daf inad in rafaranea 1 and 
Flguro 4.1 of thla annual . 

Equtclona oolvod ara Invaraa tangant funetioits and 
raqulrad logic la uaad which anauraa propar algn and 
angla quadrmc ara aalntainad. 

All angina ara ratumad In dagraaa. 

CALL ANCLES <VToTX,VTOTV,VTOTZ,VAnAL,mDAL, VTANC, 

FSI, R0TANG. 0CTA:{G,UPWASH, SIDNSH) 


Coamon Blocks 
Raqulredi 


Nota: All valocitlus given are values nomalized by 
free streaa valoclty aagnltuda. 

Inputs ; 

VkXLAL Axial velocity component In propeller plane 
coordinate system at point. 


VRADAL Radial velocity component in propeller plane 
coordinate system at point. 

VTANC Tangential velocity component in propeller 
plana coordinate system at point. 


VTOTX, X. yt and z components of velocity at the 
VTOTY, propeller plane point relative to reference 

VTOTZ body coordinate system. 

PSl Azimuth position angle of propeller plane point. 


i 

t 




' OUTaNG Outflow angle . 

ROTMG Rotational flow angle. 
UPW^H Angle of upwash. 

SIDWSH Angle of sldewash. 

Subroutines^ 

Called; None 

Error Returns: None 

Error Stops: None 








S vj) ro ut the COfeFXC 



Purpose:-/- ■' /'--To- ;C«lculate'-the\aftJ:;riX' bx'.coe^^^ the full ;~ 

■,.-"i-, ■ ..- -size-SySteatof liaeer/^equutioMV-;^^^ equations' ere/;. . ;:;-V-^ ' 

leter Wived for f he uhknot^ body'; 

Method: Ihis Subroutine Is used only .if the symetrlc inputv/ /^ 

V option has not >been --tMed. "'/"V; ■' 

An equation is generated for .eve^ ,paa«i ivvj Che 
configuration. V 


the equations (rows of natrixl are oreatW Wd d*W>dred 
sequentially starting with panel 1 then panel 2 , etc i. 


To generate an equation, the left side coefficients are 
obtained first. Then the right side coefficients are 
obtained. See reference 1 for description of the ■ 

equations. The equation is then written on one idgical" 
record of the sequential scratch file 9. 


The subroutine returns with one equation written for 
each panel and stored sequentially on file 9. 

Use: CALL COEFIC(LLL,MD,NP,XC,YC,ZC,S,ANVX,ANVY,ANVZ,A) 

Common Blocks 

Required: INPUTS, INLET, and WINGl 


Input ; 

LLL Program dimension size variable specified in 

main program which equals the maximum allowed 
number of body panels which may be generated. 
Used here for execution time dimensioning of 
arrays. 

MD Program dimension size variable sepcified in 

main program which is sum of LLL and maximum 
allowed number of body orientations. Used here 
for execution time dimensioning of arrays. 

NP Total number of body panels on configuration. 

Same as the number of equations generated. 

XC,YC, Arrays of coordinates of each body panel control 

Zt point relative to reference body coordinate 

system. 

S Array of body panel surface areas. 
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Arrays of conpoaents of the outwird drawn^^^ 

QonBBl unit vectors located 
control point* Helative to reference body 
coo.fdinate eystea._ ■ ' 

^IKFLO Total nu]d>er of body paneli specified as being 
inlet or outlet panels* / 

INDEX . Array of panel index nuoi)ers of body panels 
which have been specified as inlet or outlet- 
- panels . ■ 

NALFHA Total nui^er of body orientations. Sane as the 
nuodier of right side coefficients on each equa- 
tion being generated* 

NWING Wing signal variable: ■ 0 if no wing is modeled 

* 1 if wing is modeled* 

1 FLRATQ Array of specified velocity ratios assigned to 
the inlet and outlet body panels* 

Output ; (Written on file 9) 

A Buffer array containing left and right side 

coefficients of an equation. Used to write 
each equation on the scratch file 9 . 

Subroutines 

Called; WINGV (only if wing is modeled on configuration) 

Error Returns; None 
Error Stops; None 


ANVX, 

ANVY, 
ANV2 




yoiliii 
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SubrQUtine COFSYM ' " 

Purpose:" To calculate the natrix of coefficients of the half 

size system of linear equations. The equations are 
solved later for the unknown body panel source strengths* 

Ifethod: This subroutine is used only if the symmetric input 

option has been used.* (Body symmetry exists.) 

An equation is generated only for the panels on the left 
side of the body. The right side image panels are auto- 
matically accounted for. 

The equations (rows of matrix) are created and nun^ered 
sequentially starting with first left side panel then 
the second left side panel, etc. 

To generate an equation, the coefficients on the left 
side of the equation are obtained first. Then the 
right side coefficients are obtained. Ste reference 1 
for a description of these modified symiaelric body equa- 
tions, T’ne equation is then written on one logical 
record of the sequential scratch file 9. 

The subroutine returns with one equation written for each 
panel on the left side of the body. Tlie equations 
are stored sequentially on file 9. 


Use: ' CilLL COFSyM(LLL,MD,NP,XC,YC,ZC,S,ANVX,AMVy 

Common Blocks 

Required: . INPUTS, INLET, SYl-fTRE, and WINGl 

Input : 

LLL Program dimension size variable specified ir. 

main program which equals the maximum allowed 
number of body panels which may be generated. 
Used here for execution time dimensioning of 
arrays . 

MD Program dimension size variable specified in 

main program whic'n is sum of LLL and maximum 
allowed number of body orientations. Used here 
for execution time dimensioning of arrays. 


NPNSYM 


NP 

XC, VC, 
ZC 

S 

ANVX, 

ANVY, 
ANV2 

NINFLO 

INDEX 

NALPHA 

NIJING 

FLRATO 

Output ; 
INSOLV 


, Bookkeeping- wKieh dest^ates leffc side 
panels and recor la cortesponding ?lght side • 
...-p4Uiels-'-as;-folloWs : ' 

if NPNSVM(I) • P , then panei;^ whose Index 
nunber is I is on the i^ght side of the body. 

if NPNSYMCt} - (N d CO, then panel 

ntdbet is I is on the left side of the body 

and the corresponding panel on the right side 

of the body has index nunber r 

Total nuober of body panels on configuration. 
Equals twice the nuiiber of equations generated. 

Arrays of coordinates of body panel control 
points. (Only those on left side panels are 
used.) 

Array of body panel surface areas . 

Arrays of components of the outward drawn 
normal unit vector at each panel control point. 
Components relative to reference body coordinate 
system. 

Number of body panels which are specified as 
being inlet or outlet panels. 

Array of panel index nunbers of body panels 
which have been specified as inlet or outlet 
panels . 

Total number of body orientations. Same as the 
number of coefficients on the right side of 
each equation. 

Wing signal variable: * 0 if no wing is modeled 

= 1 if wing is modeled. 

Array of specified velocity ratios assigned to 
the inlet and outlet body panels. 


Bookkeeping array which stores the index 
numbers only of the panels on the left side of 
the body. Index numbers are stored in sequence 
such that INSOLV(l) stores the index of the 
first left side panel, INS0LV(2) stores the 
index number of the second left side panel, etc. 




Subroutines 
Called : 

Error Returns: 




Other Output ; (Written on file 9) 

A Buffer array containing left and right side 

coefficients of an equation. Used to write 
eacn equation on the scratch file 9. 


WINGV (only if wing is modeled on configuration) 
None 


Error Stops : 


None 
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Subroutine EHLER 
Purpose; 


‘i 


Method: 


Use: ' X 

Coiiunon Blocks 
Requited:^ 


To calculate valusi of an euler eagle used in coordinate': 
triuMfoti&Stion, The trana the 

refe rence body -coordinate syst'eo either the wind . 
axis or propeller fixed a^s ayaten^ 

'Solves theVinvezae tangent function expression for the 
euler angle as given in reference 1.; Logic is used to 
keep track of sign so proper quadrant value of angle is 
not lost. 




CAI^ mER(A,B»MADA) 

None , 

Input : 


',<ir 


Value of either a body angle of attack or ' 
angle of attack of the propeller plane relative 
to the body (Degrees) 

Is either the body sideslip angle or the 
yaw angle of the propeller plane relative to 
the body (Degrees) . ' 


Output ; 



ANADA 

Subroutines 


Called: 

None 

Error Returns : 

None 

Error Stops: 

None 


Euler angle used In coordinate transfornatlon . 


■’■"fr- 
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Subroutine INPUT 


Purpose: To read in all input data to the computer program. 

Method: Eight seta of data are read. (Refer to section 3 of 

this manual.) As each data card is read, the values 
are also printed along with the variable name. 

The quantities read are checked for detectable errors 
which will be explained below in the section on error 
stops. When any error is detected in a value, a message 
describing the error is printed. The subroutine will 
stop processing data and will cause the entire program 
to terminate. If no errors are detected in the data, 
the input will be saved In arrays and variables. Then 
a llsrlng will be printed which summarizes what further 
output tables to expect from the remainder of the program. 

Control is then returned to the main program. 

Use: CALL INPUT(LLL,MAXALP,^!AXINF,NIP^^AX,NSECMA,NSECTO, 

NCOUNT , NFLAG , PX , PY , PZ , ITMAX , ERR) 


Common Blocks 

Required: TITLE, INPUTS, OPTION, PROP, SYMTRE, INLET, and WINGl 

Input : (Supplied from main program.) 

NSECMA Program dimension size variable specified in 
main program which equals the maximum number 
of body cross sections which may be used to 
describe the geometry. Used here for error 
checks and execution time dimensioning of 
arrays. 

NIPMAX Program dimension size variable specified in 
main program which equals maximum number of 
periphery points which may be input around any 
cross section. Used here for error checks and 
execution time dimensioning of arrays. 

MAXALP Program dimension size variable specified in 
main program which equals the maximum allowed 
number of body orientations. Used here for 
error checks and execution time dimensioning 
of arrays. 
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LLL Progran dltaetwidn « if a yar labia aH^tflad In > 

' oai^ naxlaua poaalbla 

nuabar of body panaXa vhi?b nay ba ganaratad lby- - 
tha input i Uaad hara for atror ehac^ and 
axaeution tlsa dlmanaioning of arraya. 

HAXINF Progran dliaanalon siea variabla apaclfiad in 
aain program which aquala tha maximum numbat 
of inlat and ouclat panala chat may ba pra~ 
acribed ovar tha anclra body. tJaad hara for 
error chacka and axccution tiaa dimanaioning 
of arraya."'. \ 

Remaining Input ; (Supplied from input card deck.) 

SY^OL Array of alphanumeric charactara making up Che 
input heading . 

V Magnitude of free stream velocity. 

NALPHA Number of body orientations to be run. 

ALPHA Array of body angles of attack. 

BETA Array of body side slip angles . 

NSYMET Symmetric body input option: * 

* 0 if body symmetry is used. 

■ 1 if body symmetry is not used . 

NSECTO Total number of body cross sections used to 
define body . 

NSEC Cross section sequence number. 

NIP Number of input points on the cross section. 

NEND Cross section flag signal: 

* 0 if cross section is not repeated or not the 
end of a panel network. Place panels between 
this section and Che next. 

> 1 if cross section is repeated or is the end 
of a panel network. Do not place panels between 
this section and the next one. 

X,Y,Z Coordinates of a cross section perlpher;;/ point 

relative to the reference body coordinate system. 

NLIST Panel geometry table print option; 

■ 0 if all body panel geometry Is to be printed in 
table during execution of paneling subroutine. 

* 1 if geometry is not to be printed. 


NCALC 

NINFLO 

INDEX 

FLRATO 

NWING 

XQR, 

YQR, 

ZQR 

SPAN 

CHORD 

DIKED 

SIVEEP 

a 

NPOINT 


XHUB, 

YHUB, 

ZHUB 

AI.PH? 

BETAP 

RADIUS 


CaoiMtry check— tMt run option: 

• 0 for a normal run to otbtain flow solutions « 

■ 1 for cast run to obtain only panel geometry. 

Number of panels specified as being inlet or 
outlet panels. 

Array of panel Index nt^inbers to identify body ^ 
panels. ' -' •[ 

Array of specified inlet or outlet flow velocity 
ratios.' ^ ’■ ' 

Wing modeling signal variable: 

• 0 if no wing is modeled on this configuration* 

• 1 if wing is modeled . 

Coordinates of the wing root quarter chord 
location relative to reference body coordinate 
system. 

Physical wing span. 

Wing root chord length* 

Wing dihedral angle* 

Sweep angle of wing quarter chord line. 

Array of wing lift coefficients. 

Propeller plane signal variable: 

“ 0 if no propeller plane is modeled and no* 
propeller plane flow field predictions are made. 

■ 1 if propeller plane is modeled and flow field 
predictions will be made * 

Coordinates of the hub or center of the 
propeller plane disk relative to the reference 
coordinate system. 

Angle of attack of propeller plane relative to 
body longitudinal axis . 

Yaw angle of propeller plane relative to 
body longitudinal axis . 

Reference scale radius of propeller plane. 

Azimuth angle spacing between radial rows of 
points in the propeller plane. 


DPSI 
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Subroutines 

Called: 

Error Returns: 
Error Stops: 


DSADUS lUdlal spacing betvccn points in ths propallsr 
plana as fraction of raf trance radius. 

KRAO Nuobar of points positioned radially along the 
propeller plane. 


NPUMCH Propeller plane flow—punched card output 
option: 

■ 0 if no punched card output is desired. 

■ 1 if desire punched card output. 

ITMAX Maxiaua nuaber of iterations allowed for 
obtaining source strength solutions. 

ERR Power of ten of the value of desired solution 
error residual. 

Output ; (This output is printed only but not saved.) 


J Body orientation case number. 


Remaining Output : (All input card data is printed, but 

some of above Input values are renamed as new output 
quantities given below.) 

NIPMAX Same value as defined in input section above 
unless symmetry option NSYMEX ■ 0 was used. 

If symmetry was used, this variable takes the 
value of the maximum nusber of periphery points 
allowed to be specified on the left half of any 
cross section. 

NCOUNT Array whose 1 th value stores the NIP value 
associated with the 1 th input cross section. 

NFLAG Array whose 1 th value stores the value of 
NEND associated with the 1 th Input cross 
section. 

PX, PY, Two dimensional arrays containing coordinates 
PZ of all points (X, Y, Z) (defined In Input above). 

First array index corresponds to the cross sec- 
tion number. Second array index corresponds to 
the order in which points are placed around the 
cross section. 


None 

None 

Forty-two detectable errors will cause run termination. 
In general these errors occur if a value is out of its 
pertr.issible range. Other errors can be detected if some 
ge'' .lie try input rules have been violated. 


The dcctetabl* error* are giver below: 
Stop Number Error Dcfcriptlon 


1 V 5 0 

2 NALP2IA > MAXALP 

3 KALPHA < 1 

A NSYMET jt 0 and f> 1 is invalid. 

5 A BETA angle la nonzero but aysmetry option 
has been specified. 

6 NSECTO <3 

7 NSECTO > NSECMA 

8 NSEC not In increasing sequence starting with 1. 

9 NIP < 3 with symmetry option chosen . 

10 NIP < A with symmetry option not chosen. 

11 NIP > NIPHAX for this section— -with symmetry 

opcicn used. 

12 NIP NIPMAX for this section—with symmetry 
option not used . 

13 NEID - NFLAG(l) - 1 for first section Is invalid. 

lA NEND ft 0 on the first section is invalid. 

15 First input point on the cross section is not 
repeated as last point of Si^ctlon. This is 
required for non symmetry option. 

16 SysBsetry option used but first point on section 
is not on upper centerline of contour. 

17 Symm>^try option chosen but last point on section 
Is n?t on lower centerline cf contour. 

18 r>ynmetry option chosen but a point on the right 
side of the body has been input . 

19 NEND 0 and i is Invalid. 

20 NEND (i.e. NFLAG) « 1 specified on adjacent 
cross sections—invalid. 

21 On two cross sections between which panels are 
to be generated, the number of points on one 
section is not the same as number on the ocher 
section-invalid. 

22 Error in two sections between which panels are 
to be generated. The cross sections intersect 
or overlap since a point on one section has been 
also placed on the other. 

23 NLIST 1 and f* C is invalid. 

24 NCALC 1 and 0 invalid. 

25 Test run specified (NCALC » 1) but inflow panels 
have been specified (N’'NFL0 f 0)— invalid . 

26 NINfLO is negative. 

27 Sytametric option chosen but NINFLO value is 
incorrect (MINFLO > MAXINF/2 was used.) 

28 Non-symmetry option chosen but. NINFLO value is 
incorrect (NINFLO > JllXINF was used.) 
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Warning 
Messages : 


29 

INDEX < 1 or > LLL was used and Is 

Invalid. 

30 

NWING ^ 0 and 1 Is invalid. 


31 

Synmetrlc option was used but wing root Is not 


located on the plane of symawtry. 
YQR ft 0 was given and is invalid.) 

(i.e. 

32 

SPAN < 0 is invalid. 


33 

CHORD d 0 Is invalid. 


34 

NPOINT 0 and 1 is invalid* 


35 

RADIUS < 0 Is Invalid. 


36 

DPSX ^ 0 Is Invalid. 


37 

DPSI > 360.0 is Invalid* 


38 

DRADUS < 0 Is Invalid. 


39 

NRAD < 0 Is Invalid. 


40 

NRAD > NRADioax Invalid* 


41 

NPUNCH 0 and 7 * 1 Is invalid. 


42 

ZTMAX ^ 0 is invalid . 



Two possible warning messages, If printed, do not 
terminate program but Indicate possible problem: 

1. Warning due to DRADUS ■ 0 

2. Warning due to ERR > 0 which gives solution error 
residual greater than 1.0 
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S<rt>toutin« PASEL 


?vxpo»*t 7& f«d«rat« ch« body surface paaels from cho Input 

cross section geometry. Information produced includes 
panel Index numbers, surface areas, outward normal unit 
vector coq^onents, control points, and a total of the 
number of panels generated. 

Method: Rings of Panels are produced around the surface between 

adjacent body cross sections. The first ring is placed 
between input sections one and two, etc. The panels 
are assigned consecutive identifying index numbers 
starting with the first panel on the first ring and 
increasing back on the configuration. Within each ring 
of panels, the top left side panel is produced first 
and has the smallest index number found on that ring. 

Then panels are generated and nuabeted in sequence on 
Che ring around the left side of the body, and up the 
right side with the last pond of the ring being on the 
Cop right jide. 

Panels generated may be either triangular or four-cornered. 
A triangular panel is made when an input cross section 
point has been repeated. Otherwise, the pend will have 
four different corner points. The subroutine checks 
chat a panel of zero area has not been input. Tliis 
would happen if correspvf^ng points on both cross 
sections were both repeat^- Such errors cause complete 
program termination. The logic used by this sub routine 
for trianguJai’ panels differs from that used with four- 
cornered panels. 

With triangular panels, the input corner points arc 
coplanar an'l need not be shifted to produce a flat panel. 
The outward rortsal unit vector for this panel is obtained 
by takin,4 the vector product of the two adjacent sides 
of Che triangle. Then surface area is computed using 
the formula for the area of a triangle. Lastly, the 
control point is located at the area centroid of the 
triangular panel. 

Logic required with fovir-comereu panels Is more involved 
since the input comer points are usually not coplanur- 
First the uric normal vector is found by taking the 
vector product of the diagonals between onposite comer 
points. Nov Che input comer points are shifted to make 
them coplanar. This is done by projecting the input 
points onto a plane which contains the average of the 
input comer points and which is orthogonal to cht unit 
normal vector. The new projected points are t.hen coplanar 
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and define a flat quardilateral panel. Now the surface 
area is obtained using using the area formula for any 
general plane quadrilateral. Lastly, the control point 
is placed at the area centroid of the quadrilateral panel 

If the s>Tiiiiietric body input option, NSYyZT = 0. is used, 
only the left side of each cross section ■'•npat. 

However, the subroutine computes the ri/ht ~>iue bcd> 
panels automatically using symmetry and folicv’-: the same 
panel numbering sequence as would exist i' the geometry” 
had been input without using the symmetry option. 

For the symmetric case only, '* bookkeeping system array 
is generated which identifies a panel as bt.ing on the 
left side or on the right side of the body. This book- 
keeping array is not used if the non-symmetric input 
option, NSYMET = 1, has been used. 

At the end of subroutine execution, the total number of 
panels generated, NP, is stored and printed out. 

Finally, thio subroutine w”ill print a detailed panel 
geometry* listing if requested by the, NLIST = 0, param- 
eter. '^his listing is automatically printed if the test 
run option, NCALC = 1, has been input. This listing 
identifies all panels of the configuration. Tae panel 
index number, surface area, control point, and input 
corner points are listed. Additionally, for symmetric 
configurations, a message denotes panels which are on 
the right side of the body and states the index number 
of the corresponding left side panels. 

Use : CALL PANEL (NI PMAX, NSEOLA , NSECTO , NFLAG , NP , NCOUNT , LLL , PX , 

PY , PZ , STOXl , STOYl , STOZl , ST0X2 , ST0Y2 , ST0Z2 , S , ANVX , AN\T , 
AMVZ,XC,YC,ZC) 


Common Blocks 

Required: SY:ITRE and OPTION 

Input : (All quantities are relative to the reference 

body coordinate system.) 

NSECMA Program dimension size variable specified in 
main program which is the maximum number of 
body cross sections which may be used to 
define the configuration. Used here for execu- 
tion time dimensioning of arrays. 

NIPMAX Program dimension size variable specified in 
main program which equals the maximun number 
of points that may be used to define any body 
cross section. Used here for execution time 
dimensioning of arrays. 







/ 




V- / ■ 
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define 


>■ ?■• 


a flat quatdilateral panel . Now the surface 
^^ ^rea -ia-i^altted using using the area formula for any 
^^^^ 4 ^^^neral plane quadrilateral. Lastly, the control point 

placed at the area centroid of the quadrilateral panel 


f the s]rmmetrlc body input option, NSY>!£T » 0, is tised^ 
psnly the left side of ea<^ cross, section was . -'‘npuc. 
Owever, the sid>routine conq>utes the rii'iht ‘^loe body 
anels automatically using symmetry and follows che same 
Ipanel nusbering sequence as would exist i' the geqcetry 
had been input without using the symmetry option. 

For the symmetric case only, r. bookkeeping system array 
is generated which identifies a panel as being on the 
left side or on the right side of the body. This book- 
keeping array is not used if :the non-s}nnmetric input 
option, NSniO 1, has been used. 

At the end of subroutine execution, the total number of 
panels generated, NP, is stored and printed out. 







Use: 


/ 


Finally, thio subroutine will print a detailed panel 
geometry listing if requested by the, NLIST =0, param- 
eter. This listing is automatically printed if the test 
rvin option, NCALC - 1, has been input. This listing 
identifies all panels of the configuration. The panel 
index nunber, surface area, control 'point, and input 
comer points are listed. Additionally, for symmetric 
configurations, a message denotes panels which are on 
the right side of the body and states the index nunber 
of the corresponding left side panels. 

CALL PANEL (NIP>L4X, NSECMA,SSECTO , NFLAG, NP ,NC0UNT ,LLL , PX, 
PY , PZ , STOXl , STOYl , STOZl , ST0X2 , ST0Y2 , ST0Z2 , S , ANVX , ANVY , 
ANVZ,XC,YC,ZC) 


Common Blocks 

Required: SYMTRE and OPTION 

Input ; (All quantities are relative to the reference 

body coordinate system.) 

NSECMA Program dimension size variable specified in 
main program which is the maximum number of 
body cross sections which may be used to 
define the configuration. Used here for execu- 
tion time dimensioning of arrays. 

NIPMAX Program dimension size variable specified in 
main program which equals the maximum number 
of points chat may be used to define any body 
cross section. Used here for execution time 
dimensioning of arrays. 


91 


LLL Program dimension si^e variable specified in 

main program which equals the maximum number 
of panels which may possibly be generated on 
any conf iguration. Used here for execution 
time dimensioning of arrays. 

XSY>r£T Symmetric input option variable: 

= 0 if body symmetric input option used. 

= 1 if body symmetric input not used. 

XCALC Test run option: 

= 0 if run is normal run for flow calculations. 

= 1 if test run is made for panel geometry only. 

NLIST Panel geometry print option: 

= 0 if desire printing of geometry table. 

= 1 if do not want geometry table printed. 

NSECTO Total number of body cross sections input on 
this configuration. 

M FLAG A r r ay of input signal va r i ab 1 c s as s i on e d to 

each input cross section. I th value is 
associated with the I th cross section: 

= 0 if cross Section is not end of bodv panel 
network, or if this section is not repeated, 

= 1 it cross section is end of network, or if 
this section is repeated. 

b’COl*NT Array whose I th v.ilue equals total number of 

input points specified on the I c!i cross section. 
This is array of .\I? values for each section. 


STGXl , 

Intermediate st 

SiOYl, 

nates of input 

STOZl 

sections he twee 


being generat ed 


Intermediate st 

sTov:, 

nates of input 

ST0Z2 

se c t i on s be twee 


or two cross 


iate storage arrays containing coordi- 


o: tv:o cross 
oanels is being 


venerated. 


px, rv. 


Two dinensicnal arrays containing originw^l input 
cross section points. First arrav index equals 
cross section number of point, second index 
g-vos the sequence number of the point around 
its cross section. 














J- 









4 


» 
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Output ; (All coordinates and components relative to 
reference body coordinate system.) 

NPNSYM Bookkeeping array established only if NSYMET = 0 
to record right side panels: 

if NFSY'M(I) = 0, body panel number I is an 
image panel on the right side of the symmetric 
body. 

if NFNSi -vD = (N ^ 0), body panel 1 is on the 
left side of the body and its right side Image 
panel is number N. 

ANVX, Arrays of components of the outward dra^^m unit 
AN\^r, normal vector of each panel. I th value is 
ANVZ for I th panel. 

S Array of body panel surface areas. I th value 

is area of I th panel. 

XC, YC, Arrays of coordinates of body panel control 
ZC points. I th value is for I th panel. 

I Do loop variable which equals index number of 

panel currently under consideration. Used only 
for printing geometry. Value not saved or 
passed to remainder or program. 


STOXl , 
STOYl, 
STOZl, 
ST0X2, 
ST0Y2 , 
STOZ2 


Same as described in input. These original 
input corner points are printed on the output 
geometry table. The values are not saved after 
printing or after subroutine has completed 
execution. 


PRXl, Temporary printing output variables containing 
PRYl, coordinates of the first corner point of a right 
PRZl side image panel of a symmetric conf igurat ion . 

Xot saved after printing on geometry table. 


PRX2, 

PRY2. 

PR22 


Same as above but second corner point of 
panel. 


PRX3, 
PRY 3, 
PRZ3 


Same as above but third corner point of panel. 


c - ^ 
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Output ; (All coordinates and components relative to 
reference body coordinate system.) 

NPNSYM Bookkeeping array established only if NSYMET « 0 
to record right side panels: 

if NPSYM(I) » 0, body panel number I is an 
image panel on the right side of the symmetric 
body. 

if NFNSi -.1) » (N 5 * 0) , body panel 1 is on the 
left side of the body and its right side image 
panel is number N. 

ANVX, Arrays of components of the outward dra;<m unit 

ANVY, normal vector of each panel. 1 th value is 
ANVZ for I th panel. 

S Array of body panel surface areas. I th value 

is area of I th panel. 

XC, Y’C, Arrays of coordinates of body panel control 
ZC points. I th value is for I th panel. 

I Do loop variable which equals index number of 

panel currently under consideration. Used only 
for printing geometry. Value not saved or 
passed to remainder or program. 

STOXl, Same as described in input. These original 
STOYl, input corner points are printed on the output 
STOZl, geometry' table. The values are not saved after 
ST0X2, printing or after subroutine has completed 
ST0Y2, execution. 

STOZ2 

PRXl, Temporary printing output variables containing 

PRYl, coordinates of the first corner point of a right 
PRZl side image panel of a symmetric configuration. 

Not saved after printing on geometry table. 

PRX2, Same as above but second corner point of 

PRY2, panel. 

PRZ2 

PRX3, Same as above but third comer point of panel. 
PRY 3, 

PRZ3 


C - ^ 
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rlLw, Sanits as above but the fourth corner point of 
PRY4, panel. 

PRZ-4 


Subroutines 
Called : 


Total nurnber of body panels generated bv sub- 
routine from user input cross section geometrv. 

PX, PY, Two dimensional arrays containing all corner 
points of panels. Values in these output 
arrays are the modified values of input corner 
points if points are on quadrilateral panels. 
(Note t::ese Vvilues are not original input corner 
points and are not printed on geometry tcible.) 


STOPER Error signal 

= 0.0 if no fatal paneling errors found. 

= 1.0 if no one or more fatal paneling errors 
were found. See error stops below. 


None 


Lrror Returns: None 


Error Scops: 


If fatal error is detected, STOPER is set =1.0. Sub- 
routine continues generating panels. After alJ panel 
geometry has been handled, and if STOPER =1.0, the 
suoroutine prints fatal error m;essage and number of 
panels . wub routine aborts entire pre^gram. 


Tvpes of fatal errors detected are one of tiiree: 

1. Zero area panels have been generated due to repeat 
of bo til CO rrospond lag points on adjacent cross 
sect ions . 

d. Roundoff error in a ver.* small panel has led to 

negative panel area or t)Ciier impossible numerical 
val ue . 


3 . 


A tier mo d i f 1 c a t i on o f 
points, the result ing 
so close as to produce 
or imaginary sur ‘ace a 
a X s o • 


four- :orne red panel corner 
corner points are on a line or 
a quadrilateral having corn 
r e a . Tl: is is a r o un do f f error 
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PRX4. Same as above but the fourth corner point of 
PRY4, panel. 

PRZ4 


NP Total number of body panels generated by sub- 

routine from user input cross section geometry. 

PX, PY, Two dimensional arrays containing all corner 

PZ points of panels. Values in these output 

arrays are the modified values of input corner 
points if points are on quadrilateral panels. 
(Note these values are not original input corner 
points and are not printed on geometry table.) 

STOPER Error signal 

= 0.0 if no fatal paneling errors found. 

= 1.0 if no one or more fatal paneling errors 
were found. See error stops below. 

Subroutines 
Called: None 


Error Returns: None 


Error Stops: If fatal error is detected, STOPER is set = 1,0. Sub- 

routine continues generating panels. After all panel 
geometry has been handled, and if STOPER = i.O, the 
subroutine prints fatal error message and number of 
panels. Subroutine aborts entire program. 

Types of fatal errors detected are one of three: 

1. Zero area panels have been generated due to repeat 
of both corresponding points on adjacent cross 
sections . 

2. Roundoff error in a very’’ small panel has led to 
negative panel area or other impossible numerical 
value. 

3. After modification of four-cornered panel corner 
points, the resulting corner points are on a line or 
so close as to produce a quadrilateral having zero 
or imaginary surface area. This is a roundoff error 
also. 
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S ub rout Ino SOLSV M 

Purpose: To solve the half size system of simultaneous linear 

equations for the unknown bodv panel source strengths. 
System was generated using left side panels as the body 
geometry has been input as a symmetric body. 

Method: This subroutine is called only^ if the symmetric body 

input option was used. Tlie method of Causs-Seidel itera- 
tion is used. The system of equations is read from 
auxiliary file 9. 

The unknown source strengths on the left side panels are 
solved, and sets of unknowns for each body orientation 
are obtained. During the solution process, all input 
body orientations are being solved together during oiu^ 
execution of the subroutine. 

SoU.tion progress is printed at the end nZ each iteration. 
Residual solution error is indicated and a message 
indicates if the solution set for a body orientation has 
converged. 


Iterations continue until convergence is achieved for 
all separate sets of source strength unknowns. Should 
one or more sets of unknouTis fail to converge within the 
allotted amount of iterations, the subroutine stops the 
program and all solutions are lost. 

After the source strengths have been solved for the left 
side body panels, the proper value of source strength is 
assigned to each mirror image panel on the right side of 
the body. This is done using a panel bookkeeping array 
which keeps track of which panels are on the left and 
which are on the right side of the body. 

The solution array contains sets of source strengths, 
one set for each body orientation. Each set has one 
source strength for each panel on the body. 


CALL SOLS LLL , LLLILAF , ^ID , >L\:aL? , NT , A, 3 IC , VARIA5 
ERR , imX , S I G>L\ , S r GS AV ) 


.Sl'M, 


Co mmi'ii Slocks 
Reqtii red: 


INPUTS and SYTHRE 


Subroutine SOLSYM 


Purpose ; 
Method: 


Use: 

Common Block 
Required: 


To solve the half size system of simultaneous linear 
equations for the unknown body panel source strengths. 
System was generated using left side panels as the body 
geometry has been Input as a synsetric body. 

This s«i} routine is called only if the symmetric body 
input option was used. The method of Gauss-Seidel itera- 
tion is used. The system of eq\iations is read from 
axixillary file 9. 


The unknown source strengths on the left side panels are 
solved, and sets of unknowns for each body orientation 
are obtained. During the solution process, all input 
body orientations are being solved together during one 
''execution of the subroutine. 


_-Si»^tion progress ik printed at the end each iteration 
Hesidual solution error is indicated and a messag^—'^^^^S^ 
indicates if the solution set for a body orientation has 
converged . 


Iterations continue until convergence is achieved for 
all'Separate sets of source strength unknowns. Should 
one or more sets of unknowns fail to converge within the 
allotted amount of iterations, the subroutine stops the 
program and all solutions are lost. 


After the source strengths have been solved for the left 
side body panels, the proper value of source strength is 
assigned to each mirror image panel on the right side of 
the body. This is done using a panel bookkeeping array 
which keeps track of which panels are on the left and 
which are on the right side of the body. 


The solution array contains sets of source strengths, 
one set for each body orientation. Each set has one 
source strength for each panel on the body. 


CALL SOLSYM(LLL,LLLH.AF,MD,>LA.X.\LP,NP,A,BIG,V.\RIAB ,SUM, 
ERR, mL\X, SIGMA, SIGSAV) 


INPUTS and SYMIRE 
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Inpu t : 

LLL Proj^ram dimension specified in main program 

which is maxinum allowed numbeiT of bodv panels. 
Used here for execution time d.' monsioning of 
arrays . 

LLLUAF Program dimension specified in main program 
equal to h;i1£ the maxinun allowed number of 
body panels. Us jd here for execution time 
dimensioning of arrays. 

M‘V!C\LP Program dimension specified in main program 
equal to maximum number of body orientations 
allowed. Used for execution time dimensioning 
of arrays. 

I'f'* Program dimension specified in main program 

which is sum of LLL and M\X.\Lr. Used for 
execution time dimensioning of arra>s. 

NT Total numoer of body surface panels generated 

over the entire configuration. Is twice t!\e 
number cf equations being solved. 

F.RR Power of ten of the desired solution error. 

irMVX 'Lixinum iterations allowed in which to obtain 
solution convergence. 

N.XLriiA Number of bodv orientations (sets of unknown:^ . 

dUM Array of interm:edia te sums used in solution 

process . 

A buffer arrav containing loft and right side 

coefficients of one equation in the system. 
Contents taken from one logical record of 
file 9. 

VARIAB Vector containing in te rmod: .ice source strength 
solution value from e icli set of i::'.knv"*vv ps . Used 
to test lor solution converconee. 


1 XS 01. V A r ray c on t .i i a i ii g t a o i nd o x n umb e rs o i i * h e b e d v 

panels on the left side of t!ie cor.f igurnt Lon 
as ge ne rn t ed b v c h e nanc 1 1 n v; s ub r o u t i ne . 


Input t 


LLL 


LLLHAF 


MAXALP 



NP 

ERR 

ITMAX 

NALEHA 

aUM 

A 

VARIAB 

IXSOLV 


Program din«nslon specified in main program 
which is maximum allowed number of body panels. 
Used here for e^cution time dimensioning of 
arrays . 

Program dimension specified in main program 
equal to half the maximum allowed number of 
body panels. Usjd here for execution time 
'-dimensioning of arrays. 

Program dimension specified in main program 
equal to maximum number of body orientations 
allowed. Used for execution time dimensioning 
of arrays. 

Program dimension specified in main program 
which is sum of LLL and MAXALP. Used for 
execution time dimensioning of arrays. 

Total nu4®^^^f body surface panels generated 
over the entire^onf iguration. Is twice the 
number of ^equations being solved. 

Power of ten of the desired solution error. 

Maximum iterations allowed in which to obtain 
solution convergence. 

Number of body orientations (sets of unknowns). 

Array of intermediate sums used in solution 
process. 

Buffer array containing left and right side 
coefficients of one equation in the system. 
Contents taken from one logical record of 
file 9. 

Vector containing intermediate source strength 
solution value from each set of unknowns. Used 
to test for solution convergence. 

Array containing the index numbers oi the body 
panels on the left side of the configuration 
as generated by the paneling subroutine. 
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NPNSYM Bookkeeping array which matches index number of 
a left side body panel with the index number of 
its corresponding right side panel: 

if NPNSYM(I) = 0, the body panel whose index 
number is I is a right side panel, 

if NPNSYM(I) = (N ^ 0), the panel whose index is 
I is a left side panel and its right side image 
panel is panel number N. 

Output : (This output is also printed.) 


ITMJlX 

ERR 


Same as quantity defined in input section above. 

Specified maximum permitted solution residual 
error : 

ERROR = 10.0^^^ 


ITER Number of iteration being done. 

K Body orientation case number corresponding to 

a set of source strength solutions being found 
(This variable is not saved after printing.) 

ALPHA Array of input body angles of attack, each value 
corresponding to an orientation case for which 
a set of solutions is being found. 

BETA Array of input body sideslip angles, each 

value corresponding to an orientation case for 
which a set of solutions is being found. 


O ther Output : (This output not printed by this subroutine.) 


BIG Vector containing the maximum error residual 

obtained tor each set of unknowns during one 
iteration. 


SIGSAV 


sig:i\ 


Bookkeeping vector containing only one set of 
solutions for left side panels only. Temporary 
storage for the unknowns until source strengths 
are assigned to the right side panels. 

Array containing final N/vLPHA sets of source 
strength solution*^ for all M? body panels in 
the conf igjra t ion . 


NSICNL Convergence mint signal: 

= 0 if ail solution sets have converged. 

= 1 if one or more sets of solutions have 
raileJ to converge this i titration (See helowl. 
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NPNSYH Bookkeeping array which matches index nuidjer of 
a left side body panel with the index number of 
its corresponding right side panel: 

if NPNSYM(I) * 0, the body panel whose index 
number is I is a right side panel. 

if NPNSYM(I) = (N 0)» the panel whose index is 
I is a left side panel and its right side image 
panel is panel number N. 

Output : (This output is also printed.) 


ITMAX 


Same as quantity defined in input section above. 


ERK. >; Specified maximum permitted solution residual 
error : 

ERROR = 10.0^^^ 

ITER Number of iteration being done. 

K Body orientation case number corresponding to 

a set of source strength solutions being found 
(This variable is not saved after printing.) 

ALPHA Array of input body angles of attack, each value 
corresponding to an orientation case for which 
a set of solutions is being found. 

BETA Array of input body sideslip angles, each 

value corresponding to an orientation case for 
which a set of solutions is being found. 

O ther Output ; (This output not printed by this subroutine.) 

BIG Vector containing the maximum error residual 

obtained for each set of unknowns during one 
iteration. 

SIGSAV Bookkeeping vector containing only one set of 

solutions for left side panels only. Temporary 
storage for the unknowns until source strengths 
are assigned to the right side panels. 

SIGMA Array containing final NALPHA sets of source 
strength solutions for all NP body panels in 
the configuration. 


NSIGNL Convergence piint signal: 

= 0 if all solution sets have converged. 

= 1 if one or more sets of solutions have 
failed to converge this iteration (See below) 
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Subroutines 
Called: None 

Error ?vOturus: Mone 

Error Stops: If one or more solution sets fail to converge after 

UMAX iterations (i.e. NSTGML = I), print solution 
failure message. Subroutine stops entire program. 




1 

j 

I 

1 


1 




I 


L 


* 
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Subro utine SCI.VK 

Purpoi^e: To solve Che full si::e syston of linear equations for 

the unknown body panel source strengths. System was 
generated using all body panels because body geometry 
was input without using the symmetr>' option. 

Method: This subroutine is called only if the body has been 

input as a nonsymmet ric configurnt ion . ITie method of 
Gauss-Seidel iteration is used. ITie system of equations 
is read from auxiliary file 9. 

The unknown source strengths for all panels are solved, 
and sets of unknounis for each body orientation are 
obtained. Baring the solution process, all input body 
orientations are being solved together during one 
execution of the subroutine. 

Solution progress is printed at the end of each itera- 
tion. Residual solution error is indicated and a message 
Indicates if the solution set for a body orientation has 
converged. 

Iterations continue until convergence is achieved for all 
separate sets of source strength unknowTis . Should one or 
more sets of unknowns fall to converge within the allotted 
amount of iterations, the subroutine stops the program and 
all solutions are lost. 

If all solutions converge, the solution array is returned. 
It contains sets of source s t rengt'ns , one set for each 
body orientation. Each set has one source strength for 
each panel on the body. 

Ise : CALL SOLVE(LLL , >D, MVXALP , N? , A, BIG , VARIAB, SUM, ERR, 

rTM\A,SIGM\) 

Common Blocks 
R equi rcvi : I!iPUTS 

I nnut : 

LLL Program dimension varia’nie specified in main 

p r o g r a m wh i c :i e c ua 1 s max 1 mum all ow e d n i imh or of 
b o dy p an els. Us e J hero fore xc c u 1 1 on. 1 1 m.e 
•dinens ioning of arravs. 












Subroutine SOLVE 
Purpose: 


>Iethod: 


The unknown source strengths for all panels are solved, 
and sets of unteowhs for each body orientation are 
obtained. During the solution process, all input body 
orientations are being, solved together during one 
execution of the svLbroutine. 

Solution progress is printed at the end of each itera- 
tion. Residual solution error is indicated and a message 
indicates If the Sblut^i^ set for a body orientation has 
converged. 

Iterations continue until convergence is achieved for all 
separate sets of source strength unknowns. Should one or 
more sets of unkno:^ns fail to converge within the allotted 
amount of iterations, the subroutine stops the program and 
all solutions are lost. 

If all solutions converge, the solution array is returned. 
It contains sets of source strengths, one set for each 
body orientation. Each set has one source strength for 
each panel on the body. 

Use; CALL SOLVE(LLL,MD,MAXALP,NP,A,BIG,VARIAB,SUM,ERR, 

UMAX, SIGMA) 

Common Blocks 
Required: INPUTS 

Input ; 

LLL Program dimension variable specified in main 

program which equals maximum allowed number of 
body panels. Used here for execution time 
dimensioning of arrays. 


To solve the full size system of linear equations for V 
the unknown body, panel source strengths. System was " 
generated using all body panels because body geometry 
was input without using the syimnetry option. 

This subroutine is called only if the body has been 
input as a nonsymmetric coafiguratioh. - The methodJ^^of^^^^^^^^^^^ 
Gauss-Seidel iteration is usecl._ The”system of aquatic 
is read from.auxiliary"file/'9.--r'T- ■ , ; ^ 
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yX\AL? 


MD 


NT 


NAL?IL\ 

ERR 


IT>L\X 


SU>1 

A 


VARIAB 


Outrut : 




ERROR 


ITER 

K 


Program dimension variable specified in main 
program which equals the maximum number of body 
orientation cases allowed. Used here for 
execution time dimensioning of arrays. 

Program dimension variable specified in main 
program which is the sum of LLL and MAXrVLP. 

Used here for erecution time dimensioning of 
arrays. 

Total number of body panels in the configura- 
tion which equals the number of equations in 
the system. 

Number of body orientation cases (sets of 
unknowns) . 

Power of ten of the desired solution error 
residual. 

Maximum iterations allowed in which to obtain 
solution convergence. 

Array of intermediate sums used in solution 
process. 

Bufftjr array containing left and right side 
coefficients of one equation in the system. 
Contents read from one logical record of 
file 9. 

Vector containing intermediate source strength 
solution values from each set of unknowT.s . 

Used to check solution convergence. 

(The following output ‘.s also printed out.) 

Same as quantity defined in input section above. 

Specified maxinun: permitted solution residual 
error : 

ERROR = 10.0^^” 

Number of iteration being done. 

Body orientation case number corresponding to 
a set of source strength solutions being found 
(This variable is not saved after printing.) 





j 


1 

j 

i 


\ 

{ 

] 


J 



In main 


aliWed In which to obtain 


“s MMUiE.r: 

TrSpc;^:: - time dimensioning of 

' : --r.;- =■•:-; -'5'. ; 3 ... ^ \ 3- 


-srTTTT^srcT* !-^-T- --i panels, in the c<»£igu«i«f- 

of ''equations’ in' 


HALPHA Kt^eirof body orientation cases (sets of. 

«ff ^±fie desl^•ed solution error 


‘ ■;■>, \ ‘ 

. V.. 7/v,: -^Uj|:. 


Array of intem^diate sums used in solution 


process. 


Buffer array con taihittg left and right side 
coefficients of one equation in the system. 
Contents read from one logical record of 

fiie- '9.' ■■■'■■'■-' ^ ' ■ ■■■■- 


VARIAB Vector containing intermediate source strength 
solution values from each set of unknowns. 

Used to check solution convergence. 


Output ; (The following output '.s also printed out.) 


ITIiAX 

ERROR 


Same as quantity defined in input section above. 


Specified raaximum permitted solution residual 
error: 

ERROR - 10.0^^ 


ITER 

K 


Number of iteration being done, 


Body orientation case number corresponding to 
a set of source strength solutions being found 
(This variable is not saved after printing.) 
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I 


Subroutines 

Called: 

Error iveturns 

Error Stops: 



ALPHA Array of input body angles of attach, each 

value corresponding to an orlenta*.ion case for 
which a set of solutions is being found. 

BETA Array of inpi’L body sideslip angles, each 
value corresponding to an orientation case 
for which a set of solutions is being found. 

BIG Vector containing the maxirnum error residual 

obtained for each se** of unknowns during one 
iteration. 

Other Output ; (The following is not printed by this 
subroutine . ) 

SIG>tA Array containing final XALPMA sets of source 
strength solutions for all NP body panels in 
the configuration. 

NSIGNL Convergence print signal: 

= 0 if all solution sets have converged 
= 1 if one or more sets of solutions have 
failed to converge this iteration (see below). 

None 
: None 

If one or more solution sets fail to converge after 
UMAX iterations (i.e. NSTGNL = 1), print solution 
failure message. Subroutine stops entire program. 


1 




1 

j 


i 

I 

i 




1 

tj 

1 


l 






aim 


.J 


■■ .'■ .-n:; ■■;:'■*■■'■■ '-V r, 

ALPHA Arrijf' bf Ifipttt'lsody AftgMjKbf.'al^^ 

- — -r-rr-— — VA^ 

tfiridt a set ef solutlota is being 

1E3CA: : A«»y of input bsty^ aiftiallp ’ angJef r 'titf^ 

valro correspondiog to <m orianeatlon o«te, ^Vo 
for which a set of solutions is being. fe^« — v-^ 


•--.Asiw:''-- 


yectof..:o'oi^A|£ri|^- 
obtai^f a; for 




j frhe following is not printed by this. 


Subroutine.) ..,.:i .. 

SIQIA Array containtogtfiaai NALPHA seta^ of sduMo- 
strength solutions for atlsf^^ pad*^Jtt. 
- ■ - ■the oohiiiuratio##:^*:-^'V"''-'"'%n'r/ 


NSIQiL Convergence print st^al: ; ;:,L.i:. 

■ 0 if all solution sets haVe' ccWverged" 

■ 1 if one or were sets of solutions have f ^ 
failed to cohVeige this iteration '(see below) r 


Subroutines 
Called: None 

Error i.etums: None 


Error Stops: If one or more solution sets fall to converge after 

ITMAX Iterations (i-e. NSIGNL ■ 1) , print solution 
failure message. Subroutine stops entire program. 
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Subroutine VCOMP 


Purpose ; 

To calculate sets of components of free stream velocity 
relative to the reference body coordinate system. Each 
set of components corresponds tc a separate body 
orientation. 

Method : 

Makes a conversion from the wind axis syste.ii (aligned 
with the free stream velocity direction) to the reference 
body cartesian coordinate system, the euler angle for 
the transformation is obtained using subroutine EULER. 
Reference 1 contains the transformation formulas. 

Use: 

CALL VCOMP 

Conunon Block 
Required: 

INPUTS 



Input : 



NALPHA 

Number of body orientations for which sets of 
velocity components are to be found. 


ANADA 

Euler angle required in axis transformation 
(obtained from subroutine EULER). 


ALPHA 

Array of body angles of attack, each value 
corresponding to a body orientation case. 


BZTA 

Array of sideslip angles, each value 
corresponding to a body orientation case. 


Output : 



VX, VA-, 

vz 

Arrays of frtre stream velocity components in 
tile reference body coordinate system. Arrays 
contain sets of components, each associated 
with one bodv orientation. 

Subrout ine 
Called: 

ELHER 


Error Returns: 

Non*. 


Error S tons : 

None 









w§m. 




r,<SS^rifti#iU»*'v r^'T^ J«?;di:i"'^"^^^rJ'^ 

:1ti^A :M. xatew6^^.>fe '<^d:;m]f^.^,*:^<^ 

■ l^ef ei^^:' 1 ■ coiH&4ii^-th€\^ forafalM-P -■-'^' 


CALL V(^HP 


CdBfflXHi Block 
Required r 


IIS’UTS 


Uaiiter of body orientations for trfilcb sets of' 

■,v-',;,v:;:,V; 5 - --„?«^locity c<^oi»nts..are. to,' be fo«^.'’' -■,. ->.:=.C;3;rs 

MABA Euler angle required in axis tr^^^j^paCioiu.VK^ 
(obtained from subroutine EBL£R^* :v 

ALPM Array of body angles of at6acfcfc,ead»- «aiue r 
corresponding to a body orientation case*^ 

BETA Array of sideslip angles, ea^ value 

corresponding to a body orientation 'ts^,.'vs.:i 

Output ; 

VX, VY, Arrays of free stream velocity components in 

VZ the reference body coordinate system. Arrays 

contain sets of components, each associated 
with one body orientation. 


Subroutine 

Called: 


EULER 


Error Returns: 


Error Stoos : 
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Subroutine VELOCI 


Purpose: To calculate and print out the body surface velocities 

and pressure coefficients. 

Method: This subroutine is called for each body orientation 

case. One body orientation is being evaluated at a time. 

The flow velocity on the surface of each body panel i? 
calculated. The control point is the representative 
location of the panel. 

At panel point I under consideration, the velocity 
components induced by panel I are first obtained. To 
these components are added all the velocities induced 
at panel point I by the sources at all other body panels 
J. Then if a wing is present, the induced velocity 
components induced at panel I by the wing are added. 
Finally, the components of the free stream velocity are 
added. The resulting sums contain the y, and z 
components of velocity on the surface of panel I at its 
control point. 

Using the velocity components above, the pressure 
coefficient at panel point I is obtained. 

The velocity and pressure results at panel I are printed 
oi4 the output table along with the panel index number 
and control point coordinates. Also are printed the 
body panel source strength solutions. 


The above process is repeated until results at every 
bodv panel have been obtained. 

The velocity components and pressure coefficients are 
not saved after they have been printed. 

See reference 1 for a description of the formulas for 
velocity induced by sources and vortices. 

Us e : CALL VELOC I ( LLL , , MAXALP , I CAS E , N P , X , AXvT , AXVZ , S , 

SIGMA, XC,YC,ZC) 

Common Blocks 

Required: IMPUTS, I:TLET, and WIXCl 


Input ; 

LLL Program dimension size variable specified in 

main program which equals the maximum allowed 
number of body panels which may be generated. 
Used here for exercucion time dimensioning of 
arrays . 


Hi 




















:’ "'^j^pa^<i j[i^tc^ i3:.paa^. I *re first obtained- Xo 

'/mb 1 l^'tlie. sottsi^es at all other body paaels 

; tj^^d' at^f^^l i-by ^a wl^ are adbted. 

- '-_adde£~ ~-X^~Tis^t^^' 5u;83 eisdtaltt' t& x* y» ^d z 
Ot^^^irata^-of vo^eity oa the sarface of paael I at its 

l%iii^>l^~^lddLty ddiipuieQts above, the pressure 
coefficient at panel potot I Is obtained* 

■ ^e vel^lty' and pressure results at panel 1 are printed 
/^"ou i^prsfitpHrti^lli'alidhs with tl» panel index n&i»T 
;--and^"ti^fofVp^n1£f-'^^fdiMtdiv^"^ are printed the 

1|he p^a^BSB is ^^eated imtil results at -every 

coefficients are 

not saved af Ctfr ttey teve* been printed . 

See reference 1 for a description of the fomailas for 
velocity iT^uced by sources and vortices. 

'« 

CALL moci (LLL, M) ,MA3CALP , ICASS , NP , ANVX, ANVY ,ANVZ, S . 
SICMA,XC,YC,ZC) 


PJ5P-ii 






Use: 


Conmon Blocks 
Required: 


INPUTS, INLET, and WINGl 
Input : 

LLL Program dimension size variable specified in 
main program which equals the maximum allowed 
nudt>er of body panels which may be generated. 
Used here for e:«cutlon tine dimensioning of 
arrays . 






MAXAL? Program dimension si'te variable specified in 
main program which equals the maxinun allowed 
number of body orientations. Used here for 
execution time dimensioning of arrays. 

^G) Program dimension size variable specified in 

main program which is sum of LLL and >LAXALP. 
Used here for execution time dimensioning of 
arrays. 


ICA5E Body orientation case being considered during 
this call of the subroutine. Value is also 
printed out. 

ALPHA Array of body angles of attack. The ICASE th 

value in the array is used by subroutine. This 
value is also printed out. 

BETA Array of body sideslip angles. The ICASE th 

value in the array is used by subroutine. This 
value is also printed out. 

V Magnitude of free stream velocity. Value is 

also printed out. 


VYy .Arrays of components of free stream velocity 
VZ relative to the reference body coordinate 

system.. Pne ICASE th value of each array is 
used and printed cut. 

hP Total number of body panels on this configura- 

tion. Tnis is also orinted out. 


XC, YC, Arrays of Ci'crdinaLes of body panel ccncrcl 

points. These values appear on the print cut 
also. Values relative to reference body 
ccordinace system. 

S Array cf body panel surface areas. 

SIGMA Array of body panel source strengths. The 

ICASE th set of t’.ie strength.s is being used. 
Tliese are the solutions of the system of equa- 
tions and tney ar=: printird out on the output 
table . 


Arrays of components of outward drown normal 
lanit vectors acting at the par. • control 
points. Conponents are relative to reference 
body coordinate svstem. 
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>IAXALP 


MD 


ICASE 


ALPHA 


BETA 


V 


VX, VY, 
VZ 


NT 


XC, YC, 
ZC 


S 

SIGMA 


AXVX, 

AMVT, 

AXVZ 


Program dimension si?:e variable specified in 
main program which equals the maximum allowed 
number of body orientations. Used here for 
execution time dimensioning of arrays. 

Program dimension size variable specified in 
main program which is sum of LLL and MAXALP. 
Used here for execution time dimensioning of 
arrays . 

Body orientation case being considered during 
this call of the subroutine. Value is also 
printed out. 

Array of body angles of attack. The ICASE th 
value in the array is used by subroutine. This 
value is also printed out. 

Array of body sideslip angles. The ICASE th 
value in the array is used by subroutine. This 
value is also printed out. 

Magnitude of free stream velocity. Value is 
also printed out. 

Arrays of components of free stream velocity 
relative to the reference body coordinate 
system. The ICASE th value of each array is 
used and printed out. 

Total number of body panels on this configura- 
tion. This is also printed out. 

Arrays of coordinates of body panel control 
points. These values appear on the print out 
also. Values relative to reference body 
coordinate system. 

Array of body panel surface areas. 

Array of body panel source strengths. The 
ICASE th set of the strengths is being used. 
These are the solutions of the system of equa- 
tions and they are printed out on the output 
table . 

Arrays of components of outward drawn normal 
jonit vectors acting at the par. I control 
points. Components are relative to reference 
body coordinate system. 








N!CING Wing signal variable 

= 0 if no wing is modeled on this configuration. 
= 1 if wing is modeled. 

Output : (In addition to the input values above which 

were printed, all the following are printed but not 
saved after printing.) 

VRX, Normalized surface velocity components at a 

VRY, panel control point. Values have been made 

VRZ nondimensional by dividing by free stream 

velocity magnitude, V. The components are 
relative to the reference body coordinate 




system. 



VR 

Total magnitude of the nondimensional velocity 
at the panel control point. 


CP 

Pressure coefficient at panel 

control point. 

Subroutine 
Called : 

WINGV 

(Called only if wing is modeled 

on conf iguration 

Error Returns: 

None 



Error Stons: 

None 
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Subroutine 

Called: 

Error Returns: 


NI\^ING Wing signal variable 

* 0 if no wing is modeled on this configuration. 
- 1 if wing is modeled. 

Output : (In addition to the input values above which 

were printed, all the following are printed but not 
saved after printing.) 

VRX, Normalized surface velocity components at a 

VRY, panel control point. Values have been made 

VR2 nondimens ional by dividing by free stream 

velocity magnitude, V. The components are 
relative to the reference body coordinate 
system. 

VR Total magnitude of the nondime ns ional velocity 

at the panel control point. 

CP Pressure ccefficient at panel control point. 

WINGV (Called only if wing is modeled on configuration.) 
None 


Error Stops: 


None 


Subroutine VORTEX 


Purpose: 

Method: 


Use: 

Common Block 
Required: 


To calculate the x, y, and z components of velocity at 
point (x»y,z) induced by a finite length vortex filament 
of constant strength. 

Given the starting and endii., points of the filament, 
in the vector sense, the Biot*Savart lav for a straight 
vortex filament is used to compute the three velocity 
components at the desired location. 

If the point lies on the filament, the velocity components 
returned are set to zero to prevent numerical difficulties 
inherent in the Biot-Savart law. 

A maximum limit on the induced velocity magnitude is 
used which equals 20 percent of the free stream velocity. 
If the velocity induced at the point, using the Biot- 
Savart lav, is more than the ?iaximum limit, the induced 
magnitude is set equal to 20 percent of the free stream 
velocity. The velocity components returned are adjusted 
accordingly. 

CALL VORTEX(X,Y,Z, XI, Y1,Z1,X2,Y2,Z2, GAMMA, WX,WY,WZ) 


INPUTS 

Note: All coordinates and velocity components below 

are relative to the reference body coordinate 
system. 


Input : 

XI, Yl, Coordinates of the starting end of the vortex 
Z1 filament. 

X2,Y2, Coordinates of the other endpoint of the vortex 

Z2 filament. (Sense of the circulation is from 

the starting end to this end point.) 

X, Y, Z Coordinates of the point at which induced 
velocity is to be calculated. 

GAMMA Strength of the vortex filament. 

Magnitude of the free stream velocity. 


V 


Subroutine VORTEX 


Ptirpose: 

To calculate the x, y, and z coaponents of velocity at 
point (x,y,z) Induced by a finite length vortex filaaant 
of constant strength. 

Method: 

Given the starting and endii.. points of the filaaent, 
in the vector sense, the Biot-Savart lav for a straight 
vortex filaaent is used to coapute the three velocity 
coiqponents at the desired location. 


If the point lies on the filaaent, the velocity coaponents 
returned are set to zero to prevent numerical difficulties 
inherent in the Biot-Savart lav. 


A maxinum llnlt on the Induced velocity magnitude is 
used vhich equals 20 percent of the free streaa velocity. 
If the velocity Induced at the point, using the Biot- 
Savart law, is more than the .vaxiaum Halt, the Induced 
magnitude is set equal to 20 percent of the free stream 
velocity. The velocity coaponents returned are adjusted 
accordingly. 

Use: 

CALL V0RTEX(X,Y,Z,X1,Y1,Z1,X2,Y2,Z2,GAMMA,WX,WY,WZ) 

Conmon Block 
Required: 

INPUTS 


Note: All coordinates and velocity coaponents below 

are relative to the reference body coordinate 
system. 


Input : 


XI, Yl, Coordinates of the starting end of the vortex 

Z1 filament. 


X2,Y2, Coordinates of the other endpoint of the vortex 

Z2 filament. (Sense of the circulation is from 

the starting end to this end point.) 


X, Y, 2 Coordinates of the point at which induced 
velocity is to be calculated* 


GAMMA Strength of the vortex filament. 

Magnitude of the free stream velocity. 


V 
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Subroutines 

Called: 

Error Returns: 
Error Stops: 


I 

I 

Output: i 

I 

V\'X, Three vortex induced velocity components at j 

VVY, the point. * 

WZ 1 

! 


None 

None 

None 
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Subroutines 

Called: 

Error Returns: 


Output : 

WX, Three vortex induced velocity components at 

WY, the point. 

W2 


None 

None 


Error Stops: 


None 


Slihroutino VPROPS 


10 / 


Purpose : 


Method: 


To calculate the flow field, velocity compcnents and 
flow angles, at all points in the propeller plane. 

Subroutine returns control inrued lately to inain program 
if no propeller plane has been input, MPOINT ^ 0. 

Othenvise, during one call of this subroutine, the 
propeller flow neld at one input bodv’ orientation is 
found. 


Pc’ \cations in the propeller plane are defined in 
.-I'^ler plane cylindrical coordinate system 
• . — *■ U ^o tne azimuthal and radial snaci.ngs specified 

by thf' input data. 


ror each point, in turn, the following is performed. 
The point's position in terns of Che reference bodv 
cartesian ceordltiaCe systen is conpoteJ. The rointh 
waCi.al posit ..on is nondimens ional iaod b .* divldin*'^ bv 
the reference propeller plane radius value obtained’ 
from common block PROP. 


Xext, the velocity components 
reference body coordinate svs 
con tribut ions induced at the 
sources are summed. Tlien, if 
velocity induced at the point 


are round in tcrm> of th.e 
ten. Tne velocity com.ponent 
point oy tne body panel 
a wing is present, the 
by t’OiC wing horseshoe 


vortex is 


Finally, tne free streari vo loci tv 


coTiponanzs are added. Tne total flow velocity ccmpcncn 
and magnitude at the point are norm.aliaed bv dividing 
b y t h. e free stream ve loci t y ma gn i t ude . 


The point Location in both coordinate 
n o rm.i 1 1 z e c Vt' 1 o c 1 1 y c o r.p on on t s in t *n e 
coo ruinate system are printed out in 
flow field table part 1 for this body 


sys terns and 
reference body 
ropeller output 
orientation. 


u..e reicrenco vcartesian velocit** components are 
trinsfermec to normaiiced axial, radial, and tangential 
components in the proneller plane cylindrico.l axis 
system. Using the vdocitv components in both cov)rdlnatt 


s y s t c ms , t h e u p a s ; i a n l: 


\ - 


1 (.10 




at tile 


poi: 


ar.- found. Also t!ie rotational flow and outfl.v. angles 
as defined in reference 1 and Figuro -:*.l of this maniiai 
are calculated. 


.ae point s d imons ion l.'ss radial ard an^muthnl oosition, 
.1 .•** vi i , t an ge n t ^ a » , an c racial ve I o c ! t v o omo on e r' t , a 
w ell a all f an y • s :\ v e 
sc r Itch, fil^ 


in 


.rtiten to one 
"o. Phis ^inr^K'S printing r.nv 
t:i.:^e later darin t t'.i is 'a:'^*on: 


: nan * 


‘ ■it i , 
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Purpose : 


To calculate the flow field, velocity components and 
flow angles, at a?.l points in the propeller plane. 


Method : 


Subroutine returns control immediately to main program 
if no propeller plane has been input, NPOINT - 0. 

Othen^.use, during one call of this subroutine, the 
propeller flow field at one input body orientation is 
found. 

Po * ’w cat ions in the propeller plane are defined in 

.j^"^ler plane cylindrical coordinate system 

g; to the azimuthal and radial spacings specified 
by the input data. 

For each point, in turn, the following is perfoirmed. 

The point's position in terms of the reference body 
cartesian coordinate system is computed. The point's 
radial position is nondimensionalized hy dividing by 
the reference propeller plane radius value obtained 
from common block PROP. 


Next, the velocity components are found in terns of the 
reference body coordinate system. Pne velocitv component 
contributions induced at the point by the body panel 
sources are summed. Tlien, if a wing is present, the 
velocity induced at the point by the wing horseshoe 
vortex is added. Finally, the free stream velocity 
components are added. The total flow velocity components 
and magnitude at the point are nomalized by dividing 
by the free sti'eam velocity magnitude. 

The point location in both coordinate systems and 
nonnalLzed velocity components in the reference hcdv 
coordinate system are printed out in propeller output 
flow field table part 1 for this body orientation. 


Next, t he reference c a r t es i an ve lo c i t y c omp on e n c s a re 
transf:?med to normalized axial, radial, and tangential 
components in the propeller plane cylindrical axis 


s y s t e m . Us in g t h o ve lo c i t y comp on o n c s i 
s y 5 1 c. ms , the up w ash an d s i d e w a s i i a n g 1 e s 
are found. Also t!ie rotational flow' and 
as defined in reference 1 and Figure 4.1 
are calculated . 


n both coordinate 
at thiO point 
out flo'A angles , 
of tills manual. 


Th e point's d 1 me ns i o n I ^ ? s s r ad i a 1 a r d a z i mu t h. a I e o s i 
a X 1 a i , t an gent i a I , a : i -d r a c i .'i I ve 1 o c i t y o o m p o n o n t s , 
w’ e 11 as a 11 f 1 o v aai g L t.' s a r e w r i 1 1 tui ^ o on r e : o r d 
s c r a t c l\ rile c d* . T : i i s e C; a 1 c s p r 1 n t i n y t ; e s e c: u a n. t 
in a second table later caring this subroutine c-xee 


ion 


t ; on . 
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Use: 

Common Blocks 
Required : 


All of the above process is repeated until all points 
in the propeller plane have been considered. 

After all points have been considered, the contents of 
file 50 is printed in the propeller plane flow field 
output table part 2. 

Lastly, this subroutine produces card ^>unched output if 
the punch option, NPUNCH = 1, was specified. This 
punched output consists of a card for each propeller 
point giving point location and the calculated axial 
and tangential velocity components at the point. See 
the output description section of this manual for more 
details on the punched output. 

The formulas for transformation of coordinate systems, 
formulas for induced velocity by sources and vortices, 
and formulas for the flow angles are ail given in 
reference 1. 

Note that all newly calculated outout quantities from 
this subroutine are printed out or punched on cards but 
are not saved after printing or punching. 

CALL VP ROP S ( LLL , MD , M-\X.\LP , I CASE , NP , AN\^ , ANIT , AN VZ , S , 
SIG^LA,XC,YC,ZC) 

TITLE, INPUTS, PROP, OPTION, WINGl, and INLET 


Inpu t : 

LLL Program dimension size variable specified in 

main program which equals the maximum number 
of panels which may possibly be generated on 
any configuration. Used here for execution 
time dimensioning of arrays. 

MAIvALP Program dimension size variable specified In 
main program which equals the maximum number 
of body orientations which may be input. Used 
here for execution time dimensioning of arrays. 

MD Program dimension size variable specified in 

main program which equals the sum of LLL and 
MAMALP. Used here for execution time dimen- 
sioning of arrays. 

Propeller plane signal variable: 

= 0 if no propeller plane is specified. 

= 1 if oroneller olane is snecified. 


NPOINT 
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Use: 

Common Blocks 
Required : 


All of the above process is repeated until all points 
in the propeller plane have been considered. 

After all points have been considered, the contents of 
file 50 is printed in the propeller plane flow field 
output table part 2. 

Lastly, this subroutine produces card ^’‘unched output if 
the punch option, NPUNCH - 1, was specified. This 
punched output consists of a card for each propeller 
point giving point location and the calculated axial 
and tangential velocity components at the point. See 
the output description section of this manual for more 
details on the punched output. 

The formulas for transformation of coordinate systems, 
formulas for induced velocity by sources and vortices, 
and formulas for the flow angles are all given in 
reference 1. 

Note that all newly calculated output quantities from 
this subroutine are printed out or punched on cards but 
are not saved after printing or punching. 

CALL \TR0P S ( LLL , MD , MAXALP , I CASE , NP , ANVX , , AN VZ , S , 

SIGMA, XC,YC,ZC) 


TITLE, INPUTS, PROP, OPTION, WINGl, and INLET 

Input : 

LLL Program dimension size variable specified in 

main program which equals the maximum number 
of panels which may possibly be generated on 
any configuration. Used here for execution 
time dimensioning of arrays. 

>L\XALP Program dimension size variable specified in 
main program which equals the maximum number 
of body orientations which may be input. Used 
here for execution time dimensioning of arrays. 

MD Program dimension size variable specified in 

main program which equals the sum of LLL and 
^L■\XALP• Used here for execution time dimen- 
sioning of arrays. 

Propeller plane signal variable: 

= 0 if no propeller plane is specified. 

= 1 if propeller plane is specified. 


NPOINT 
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ICASE Body orientation — lift case nunher of the 

orientation for vhich the subroutine is now 
getting propeller pU'ne flew field output. 

ALPHA Body angle of attack cor responding to 
orientation case ICASE. 

BETA Body sideslip angle corresponding to orienta- 
tion case ICASE. 

NVING Wing input signal variable 

= 0 if wing is not r*odeled on this configuration. 
= 1 if wing is nodeled. 

CL Array of lift coefficients of wing (if present) 

Value used is the one corresponding to orienta- 
tion case ICASE. 

V Free stream velocity magnitude 

Va, \T, Arrays of components of the free stream 

VZ velocity. The values used are the ones corre- 

sponding to orientation case ICASE. Cemponents 
are relative to the reference body coordinate 
system. 

NINFLO Number, of body panels specified as having 
inflow or outflow th^f-ough their surface. 

XHUB, Position coorcinates of center of the prcpellur 

YHUB, plane relative to reference body coordinate 

ZHL’B system. 

RA.DIUS Refc!rence propeller plane (blade) radius. 

ALPHP Pitch angle of propeller plane thrust axis 
relative to reihrcnce bod;, x axis. 

3ETAP Yaw :!ngie f prcpeller plane thrust axis 
relative to reference bed/ x axis. 

D?SI Aaimuch angle spacing of points around propeller 

plane. 

NPCNCH Punched card output o;.'tlon variable: 

= 0 if no output is to be punched on cards 
= 1 if desire punched card output. 

SlYIBOL AlphaniirriCric array containing the heading 
supplied :!n the input data. 


« V 


m 







gfeeKJSsiee3b©dv eobitdaj^feer 


BADIIB^ ^i^efere ncsT Prop&les pi^jS- Cfeia^) 

ALPHP - Pitch '^gie &f propeller^ ; 

. relative to reference body x axis; / ~ ' 

BETitf . Yaw angle of propeller plane thrust a^s - s 

relative to reference body x axis. - 

DPSI Azimuth angle spacing of points around propeller 
plane,.- 


NPUNCH Punched- card output option variable: 

' “0 if nb' output -is to"&e'punched on cards 

- 1. if desire punched card output. 

SYMBOL Alphanumeric array containing the heading 
supplied in the input data. 








110 


DRADUS Radial spacing of points in the propeller plane 
in fraction of reference propeller plane radius 

NRAD Number of points located radially along each 
azimuth position in the propeller plane. 

NP Total number of body surface panels in the 

configuration. 


XC, YC, Array of body panel control point coordinates 
ZC in terms of reference body coordinate system. 


S 

SI OLA 


Array of 

Array of 
with the 


body panel surface areas. 

body panel source strengths associated 
body orientation case, ICASE. 


ANVX, Arrays of components of the outward drawn 
AN\^’, unit normal vector at each bod/ panel in terms 

ANVZ of the reference body coordinate system. 

Output : (All output printed and/or punched. All new 

output below is lost after printing and after the 
subroutine hvis completed.) 


Old input quantities which are also p rinted and sa ved : 


ICASE 

Same 

as 

input 

de fini tion . 

ALPHA 

Same 

as 

input 

definition . 

BETA 

Same 

as 

input 

definition . 

CL 

Same 

as 

input 

dciiinit ion . 

V 

^ime 

as 

input 

de f inition . 

VX, vT, 

vz 

Sane 

as 

input 

de fini tion . 

:vHU3 , 

yhlt;, 

ZHUB 

Same- 

as 

i nput 

de f inition . 

R.ADIUS 

Sane 

as 

input 

de f ini tion . 

AL?H? 

Same 

as 

input 

definition . 

BETA? 

Same 

a s 

input 

def ini tion . 
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FLRATO Through flow velocity assigned tu first inlet 
of outlet panel (if any) specified in the 
program input. Quantity is fraction of free 
stream velocity magnitude. 

New output no t s;^vod : 

RR.\TIO Dimensionless radial position of point in 

propeller plane in the propeller plane cylin- 
drical axis system. Is ratio of radius to 
reference propeller plane radius. 

?SI Azimuth angular position of a point in propeller 

plane. In the propeller plane cylindrical axis 
system. 

POINTX, Position of a point in the propeller piano 

POINTT, relative to reference body coordinate system. 

POINTZ 

VTOTX, Normalized components of complete flow velocity 

\TOTY , at a point in the propeller plane relative to 

VTOTZ the reference body coordinate system. 

VTOT Normalized resultant magnitude of the induced 
velocitv at a point in the propeller piano. 

Furthe r ou tp ut nr in ted using file 50 : 

VAXIAL, Normalized axial, radial and tangential velocity 

VR.-\P.\L, components at a point in the propeller plane 

VT.\XG relative to the propeller plane cylindrical 
axis system. 

XOTAXc'r Rotational flow angle at a noint In propeller 
plane. 

OUTANG Outflow angle at a point in propeller plane. 

UP'\A!^r. Angle of upwash ut a point in propeller plane. 

SIDW5H Angle of sideva-:: at a poin^. in propeller plane. 

J^unche d o a pd out p ut : (If requos tod) 


SYMBOL 

Same 

as 

input de f ini r ion 

Xiun 

Same 

as 

input definition 

lunu's 

Same 

as 

inout definition 


i 


i 
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ALPHA 

BETA 

PSI 

RRATIO 

VAXIAL, 

\n:ANG 


Same as input definition. 

Same as input definition. 

Same as defined earlier in output section. 

Same as defined earlier in output section. 

Same as defined earlier in output section. 


^TSIS Total number of azimuth locations at which 
points are situated on the propeller plane. 

ATOT Propeller plane angle of attack between propeller 
plane thrust axis and free stream velocity 
direction. 

3T0T Propeller plane angle of sideslip between 
propeller plane thrust axis and free stream 
velocity direction. 

Subroutines 

Called: EULEP. WINGV (if >?.vIMG = 1), and ANGLES 


Warriing 

Messages: If 360 is not an integer multiple of DPSI , DPSI is 

printed v;ith a warning stating uneven azimuthal spacing 
of propeller plane points exists. Subroutine continues 
CO execute normally. 


Error Returns: None 


Error Scops : 


None 
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Subroutine WGEOM 


Purpose: To set up the wing horseshoe vortex geometry and strengths 

corresponding to each input body orientation-lift case. 

Method: This subroutine is called only if a wing has been 

modeled. For each input body orientation, a horseshoe 
vortex is generated. The horseshoe vortex has a span 
equal to (0.25 * ^ * physical wing span). The trailing 
vortices are shed parallel to the free stream velocity 
and are extended aft a distance of one hundred wing 
chord lengths . 

The end point coordinates of the four finite parts of 
each horseshoe are calculated relative to the reference 
body coordinate system. These four parts are the straight 
line vortex filaments consisting of the left side 
trailing vortex, the right side trailing vortex, the 
left semispan bound vortex (on the quarter chord line), 
and the right semispan bound vortex (also on the quarter 
chord line). 

For each body o rier ^ation , a different vortex strength 
is calculated as a function of the input lift coefficients. 

The output horseshoe vortex geometry and strength is 
prihted out on a cable, and the values are saved in 
common blocks ^’INC1 and WING2. 


Use : 

CALL W 

GEOM (Only if wing has been modeled) . 

Common Blocks 
Required : 

INPUTS 

, WINGI, and UING: 


Input : 



V 

Ma gni t ude of f ree stream ve loc i t y . 


SPAN 

Physical wing span. 


CHORD 

Mean or root chord of wing. 


XQR, 

Coordinates of wing root quarter chord location. 


VnR, 

Also the inboard end points of each bound sem‘- 


ZQR 

span vortex filament. Relative to the reference 
body coordinate system. 


DIHED 

Wing and bound vortex dihedral 'V'.., 'e. 













p ^ 




; y , -^ , .- . . ]^j^|M4 

SPAN Physical wing spang •^—^■■ v 
CHORD Mean or root chord' of wing. 




Coordinates of wing root ‘quarter chord location-. 
Also the^:±nboayd_P:d points of each,bb«iid/s^i#^ 
span vortex filaraent. Selative to ..the reference 
body ..coordinate aystem.. - ^-'C:;;'\:g;'.g-- :;.gV 


DIHH) Wing and bound vortex dihedral /iftgle. 










O 


SWEEP Wing quarter chord and bound vortex sweep anv^e. 

NALPHV Number cf Input body orientations . 

CL Array of wing lift coef f icients . Each va'ue 

corresponds to a different body orientation. 

ALPHA Array of body angles of attack. Each value 
corresponds to a different body orientation. 

BETA Array of body sideslip angles. Each value 
corresDonds to a different bodv orientation. 


Output : (All of the above input values as well as output, 

are printed out. All coordinates below a^e relative to 
the reference body coc’^dinate system.) 

SPAN Horseshoe vortex span. (Output v^alue replaces 

the original input value of pliy-ical wing span.) 


XBTTPL, Coordinates of the outboard enu of the left 
YBTTPL, side bound vortex filament. Also the upwind 
ZBTIrL end of the left trailing vortex filament. 

.KBTIPR, Coordinates of the outboard end of the right side 
Y 3 TIPR, bound vortex filament. Also the upwind end of 
ZBTIPR the right tniling vortex filament. 

CA>RLA Array of horse.shoc vortex strengtlis. Each value 
corresponds to a. different body orientation. 


tiDT'^u tines 
ailed : 


’“01* ..e turns 


Error Stops: 


XTR.\LL, Arrays of coordinates 'f the. dewm/ind end cf 
YTRALL, the left trailing vortex flia.ment. Each set 
ZTIL^LL of coordinates correspoTids to a different body 
orientat ion. 

XTFLVLF’, Arrays of cooriinate.i of the downwind end of 
YTRi\LR, t;ie right trailing vortex filament. Each set 
ZTiL\LR of coo rdi nsit c*s cor rt--sp onus to a different bod" 
orientat icn. 




Subr outi ne 
Purpose : 

Method : 


Use: 

Conur.on filoc 
Required : 
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InTNGV 


To calculate the velocity components Induced at a 
specified point by the wing horseshoe vortex of constant 
strength. 

The horseshoe vortex is broken into four finite length 
straight vortex filaments all of the sane strength. The 
four parts are the left side trailing fiJ ament, left 
semispan bound filament, right semispan bound filament, 
and right side trailing filament. The end points of 
the four filaments are obtained from common blocks WIMCl 
and WING2. At the desired point, the induced velocity 
components due to each of the four filaments are computed 
in turn, by the Bi’ot-Savart law. Then the component 
contributions due to all filaments arc summed to give 
the resulting velocity components due to the horseshoe 
vortex as a whole. 

CALL WINGV( ICASE , X , Y , Z , VELX, VLLY , VELZ) 


WINGl and \:VAG2 

Note: All coordinates and velocity components below arc 

relaf^ve to the reference body coordinate system. 


Incut : 


X. Y, 

Coord inac 

es of 

th 

e 

p 

oi 

nt at 

wh 

ich 

the 

vortex 

z 

Induced v 

elocit 

ie 


a 

re 

to be 

1 

ound 

• 



ICASE 

Body orie 

ntat lo 

11“ 

- 

li 

t c 

case 

(angle 

0 f 

a 

ttach , 


sideslip. 

and 

li 

c 

t 

Ov"‘ 

of ficien 

t) b 

ein 

z 



conside re 

d 











XIPaLL, 

Arrays 

coord 

in 

a 

1 


0 f the 



s t r 

t*a 

m sta.rt 

VTR.U.L, 

end of lA 

e left 

s 

ido 

t 

ra li ’ n 

rr 

vert 

ex 

fi 

lament 

ZmiL 

for the b 

i^Jy or 
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n 

La 

Li 

on cr*s 
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r:. 



:-3TIPL , 

Coor iinat 

s- S ' .) t 
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oi 

nt. ar 
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Icf 

c 

.S i if 

YSTIPL, 

c rai ling 

70 r t e 

j 

0 

in 


the ie 

[ r 

s t. m 


an 

bound 

ZBTIP!. 

vortex :i 

lame n t 

• 











tant 

strong 
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of 

r 

ho horse 

ijtiOe 

VO 
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its four 

s r r n i s: 

l.t 


ji 

ne 
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: T'3 

rts 

. 

Va 1 no 


corresponds to 

th 

e 

c 
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ent:a:i 

on 

c as 


CA 

iZ. 
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XQR, 

YQR. 

ZQR 

Coordinates of the inboard or root ends of the 
left and right semispan bound vortex filaments 


X3TIPR, 

YB^IPR, 

ZBTIPR 

CooiJinates of the point at which the right 
side trailing vortex filament joins the right 
semispan bound vortex filament. 


XTRALR, 

YTRALR, 

ZTRALR 

Array of coordinates of tho dfr^Tistrean end 
of tne right s L ie trailin:; vortex filament. 


Output : 



VELX, 

VELY, 
VELZ 

Components of velocity induced at the input 
point by the complete wing horses’ioc vortex. 

Subroutines 

Called: 

VORTEX 


Error Returns: 

None 


Error Stops : 

None 



1 

i 

I 

i 

( 

j 

j 

t 

i 

i 


Source Progran Listing 


This version of the progran is dimensioned using Che following 
values of the dinension size variables specified in the main program. 
(Quantities are explained in sections 2.2 and 2.3 of this manual.) 

NSEC:U = 60 

NIPMA.X = 45 

M4XALP = 6 

>L4XIb’F = 500 

LLL = 2596 
>© = 2602 
LLLHAF = 1293 

and : 

D?SI . = 1.0 

irin 

NRAD = 51 

raa.\ 
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SECTION J AND THE THIRD CORNER ON SECTION J+1 : PAN 7250 

(IE. PT(J+1,K) - PT(J+1,K+1) ). PAN 7300 

PAN 7350 
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UR«SQRT( ( (ST0X2 (K+1 )-ST0X2 (K) )**2 )+( (ST0Y2 (K+1 )-ST0Y2 (K) )**2 )+ PAH 9050 
l((STOZ2(K+l)-STOZ2(K))**2)) PAM 9100 

STANG-( (RS**2 )+(SU**2 )-(UR**2 ) )/(2 .*RS*S0) PAN 9150 

PAN 9200 

CHECK IF ABS(STANG) LESS THAN OR EQUAL TO I. IP NOT, NUMERICAL PAN 9250 
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IF THIS WAS THE LAST STATION OF PANELS TO BE DONE, THEN NP PAN41450 
CONTAINS THE TOTAL NUMBER OP PANELS GENERATED FOR ENTIRE CONFIG. PAN41SOO 
THEN RETURN TO MAIN PROGRAM PAN41550 
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RETURN ANG 3150 

END ANG 3200 


APPENDIX C 


Sample Case 


C.l Configuration 

To illustrate the Input and output for the computer program^ an 
exaaqile has been run. 

A simple configuration made up of 16 triangular body panels is 
used as shown in the three-view drawing of Figure C.l. The unit length 
scale is shown on the y axis in the top view of the configuration. The 
body consists of two octahedrons in tandem but not connected. This two- 
part body is used to illustrate the paneling of bodies having several 
body panel networks. The origin of the coordinate system is at the 
nose of the body, on the plane of symaetry. Panel numbers, as will be 
generated by the program, are shown on Figure C.l where panels are 
visible. 


Body Characteristics ; (length unit shown on top view of Figure C.l) 


length overall 
length front part 
length aft part 
width 
height 


4.25 units 
2.00 imlts 
2.00 units 
2.00 units 
2.00 units 


Four panels will be designated as inlet or ourlet panels. Panels 1 
and 4, on the top front of the body, are inlet panels with an inlet 
velocity ratio of 0.1. Panels 14 and 15, on the lower rear of the body, 
are outlet panels with an outlet velocity ratio of 0.1. 


A propeller plane is positioned as shown in Figure C.l. It is 
oriented perpendicular to the body longitudinal axis. Points on the 
propeller plane are shown on the front view in Figure C.l. 


Propeller Plane Characteristics ; (length unit shown on top view of 
Figure C.l) 


hub coordinates 
reference radius 
azimuth position of points 
radial number of points 
radial spacing of points 


(0.25, 0.0, 0.0) 

1.5 units 

0, 120, and 240 degrees 
11 

10 percent of reference radius 
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A reccangtilar unswept wing having no dihedral or sweep Is attached 
to the front part of the body as shown In Figure C.l. 

Wlna Characteristics ; (length unit shown on top view of Figure C.l) 

span 5.0 units 

chord 0.5 units 

location of root qiiarter chord (-1.25, 0.0, -0.5) 

wing Cj^ 0.1096 at 0 degrees body angle 

of attack 

1.0960 at 10 degrees body angle 
of attack 

The configuration will be run at two orientations. First Is 0.0 
degrees body angle of attack and sideslip, second Is 10.0 degrees body 
angle of attack and 0 degrees sideslip. 

C.2 Input Data Card Organization 

On the following two pages is a listing of the input card data 
for this sample case. The first and last lines of each page are not 
data but Indicate card column nunhers. 

The body geometry is described using six cross sections. Since the 
symmetric Input option can be used for this configuration, only the 
left side of the sections is given. Thus, only three points are 
specified on the cross section peripheries. Note, had the symmetric 
option not been used, then five points per cross section would have 
been needed. 

This will be a normal flow prediction run since the geometry is 
simple and the panels may be Identified without making a geometry test 
run. Also, the listing of body panel geometry is requested. 

Of the four inlet/outlet panels, only the two panels on the left 
side, numbers 1 and 14, are specified. This is in accordance with 
symmetric input rules. 

Refer to the following two pages for the actual input card data. 
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C.3 Printed Output Listing 

The printed output obtained for the sample case run Is reproduced 
on the remaining 15 pages of this appendix. 




Mf^G HCI>SESHC6 VOUTEX STRENGTHS AND CfC^rTAV 
ItOO^OtNATES AEtATtVC TO bOOV AXI^ SVSTfHI 




Ai 
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